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ABSTRACT 

Two quantitative tests DIFF1 and DIFF2 for measuring goodness-of-fit 
between two locally-normalized supernova spectra are presented. Locally- 
normalized spectra are obtained by dividing a spectrum by the same spectrum 
smoothed over a wavelength interval relatively large compared to line features, 
but relatively small compared to continuum features. DIFF1 essentially measures 
the mean relative difference between the line patterns of locally-normalized spec- 
tra and DIFF2 is DIFF1 minimized by a relative logarithmic wavelength shift 
between the spectra: the shift is effectively an artificial relative Doppler shift. 
Both DIFF1 and DIFF2 measure the physical similarity of line formation, and 
thus of supernovae. DIFF1 puts more weight on overall physical similarity of 
the supernovae than DIFF2 because the DIFF2 shift compensates somewhat for 
some physical distinction in the supernovae. Both tests are useful in ordering 
supernovae into empirical groupings for further analysis. We present some exam- 
ples of locally- normalized spectra for Type lib supernova SN 1993J with some 
analysis of these spectra. The UV parts of two of the SN 1993 J spectra are HST 
spectra that have not been published before. We also give an example of fitted 
locally-normalized spectra and, as an example of the utility of DIFF1 and DIFF2, 
some preliminary statistical results for hydrogen-deficient core-collapse (HDCC) 
supernova spectra. This paper makes use of and refers to material to found at 
the first author's online supernova spectrum database SUSPEND (SUpernovae 
Spectra PENDing further analysis) Jj 
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1. INTRODUCTION 

An ordinary % 2 -like test measuring the goodness-of-fit of two supernova spectra often 
fails to be consistent with what the eye sees qualitatively: good agreement to the eye can 
be poor or moderate agreement by a x 2 -like test. Since the human eye is an excellent 
pattern recognization tool, judgments based on eye comparisons informed by a specialist's 
understanding of spectrum formation are often preferred over x 2 -like tests or other formulaic 
tests in identifying similarity and/or goodness-of-fit between two spectra. Human judgment, 
of course, is qualitative and to some degree subjective. Thus, it would be advantageous to 
have tests that measure spectrum similarity more consistently with what the eye sees and 
yet be quantitative and objective. Such tests could be applied in an automated fashion and, 
one hopes, would show correlations missed by the eye. 

In this paper, we present two quantitative tests DIFF1 and DIFF2 for measuring 
goodness-of-fit between two supernova spectra. The tests are both measures of the rela- 
tive difference in line patterns. (The formulae for DIFF1 and DIFF2 are given, respectively, 
in §§ 4.1 and 4.2.) Both tests depend on what we call local normalization which largely 
removes continuum features without distorting the line features too much. Local normaliza- 
tion reduces the continuum shape to an apparent flat line of height 1 by dividing the original 
spectrum by a version of the spectrum smoothed over wavelength interval relatively large 
compared to line features, but relatively small compared to continuum features. Ordinary (or 
global) normalization is just to divide the spectrum by a constant chosen for some purpose. 
The near-equivalent of local normalization has frequently been used for synthetic spectrum 
calculations simply by inputting a flat continuum at the base of a model atmosphere or by 
dividing a synth etic spectrum by a know n continuum. In the context of synthetic supernova 



spectra see, e.g.. lDessart fe Hillierl (120051 ) for this kind of procedure. But note that local nor- 
malization is a particular procedure applied to spectra with line features and it cannot be 
exactly equivalent to using a flat continuum or dividing a spectrum by a known continuum 
in synthetic spectrum calculations. 



DIFF1 and DIFF2 evolved from the DIFF test presented by Branch et all fl2006ri l. 
DIFF1 essentially measures the mean relative difference between the line patterns of locally- 
normalized spectra and DIFF2 is DIFF1 minimized by a relative logarithmic wavelength 
shift between the spectra: the shift is effectively an artificial relative Doppler shift. From 
this description, the reader can roughly understand both tests before knowing the test for- 
mulae. Both DIFF1 and DIFF2 measure the physical similarity of line formation, and thus 
of supernovae. DIFF1 puts more weight on overall physical similarity of the supernovae than 
DIFF2 because the DIFF2 shift compensates somewhat for some physical distinction in the 
supernovae. Both tests are useful in ordering supernovae into empirical groupings for further 



- 3- 



analysis@ 



There are several reasons for local normalization (or some similar continuum elimina- 
tion technique) in comparative supernova spectrum analysis. First, the intrinsic supernova 
continuum can often be quite uncertain. Often the main reason for this is uncertainty in the 
reddening correction. Error in reddening correction can affect, for example, the B — V color 
and other colors (particularly toward the blue) by a significant fraction of a magnitude or 
more. Uncertainty in reddening is particularly a problem for core-collapse supernovae (i.e., 
Types II- P, II-L, Iln, lib, lb, Ic, and Ic hypernovae) which tend to arise in or near star- 
forming regions in their host galaxies. Foreg round Galactic redde ning can be reasonably 
con fidently corrected for usin g the resul t s of ISchlegel et al.l (119981 ) for Galactic reddening 
and ICardelli et al.l (119891 ) and lO'Donnelll (119941 ) for the reddening law, but the host galaxy 
reddening can usually only be corrected for with great uncertainty relying on interstellar 
lines in the supernova spectra or, with less uncertainty depending on cases, on spectral mod- 
eling. For example , early spectr al modeling of Type II- P supernovae can constrain reddening 
(IBaron et al.ll2003l . 120041 . 120051 ). Nevertheless, in spectral modeling one would much prefer 
to have reddening as a given rather than as a parameter to be fitted for — or not to be a 
consideration which is what local normalization helps toward. Core-collapse supernova types 
are all rather heterogeneous types, and so judging reddening by a comparison of a supernova 
to other examples of its type is not usually conclusive. Type la supernovae (SNe la) are 
usually much less reddened as they do not preferentially arise in or near star forming regions 



and the homogeneity of most of them (e.g.. lBranch et al.ll2006bl ) allows reddening correction 



by comparison in so me cases. But even SNe la can be hi ghly reddened (e.g., SN 1986G 
(jPhillips et al.lll987l )) and there are peculiar SNe la (e.g., iBranch et al.l l2006bl ) for which 
comparisons to determine reddening can be uncertain. Thus, the reddening correction for 
SNe la is often uncertain. 

In addition to reddening, the continuum can be in error because of errors in the broad- 
band calibration of the spectrum. A particular problem in achieving good calibration arises 
because supernovae are transient, time-dependent sources: this means one usually cannot 
simply replace a bad observation by a good one for a given epoch and this limits the ability to 
achieve good calibration. Another related particular problem arises from the fact that bright 
well-observed supernovae are relatively rare. Thus, many well-observed supernovae are from 



2 While this paper was in preparation, analysis tools for supernova spectra that are in some re- 
spects similar to DIFF 1 and DIFF2 with loc al normalization were de scribed by Stephane Blondin 
(jBlondin fc Tonrvl 120071 ) and Avet Harutyunyan (jHarutvunvan et al.l I2005T ) at the conference The Multi- 
coloured Landscape of Compact Objects and their Explosive Origins, Cefalu, Sicily, 2006 June 11-24, URL: 
http:/ /www. mporzio.astro.it/cefalu2006/ . 
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earlier instrumentation epochs. One would like to use the data from these supernovae for 
analysis since they are to some degree unique even though the data may have very uncertain 
calibration by modern standards. 

Reddening and calibration errors tend to affect the spectrum continuum multiplicatively: 
i.e., they change flux level by a scaling factor that varies slowly with wavelength. These 
kind of uncertainties are ideally, and nearly practically, irrelevant to the process of local 
normalization (see § 2.1), and thus cannot much affect analysis with the locally-normalized 
spectra. 

A kind of continuum uncertainty that is not eliminated by local normalization is the 
additive uncertainty due to contamination by extraneous sources. (See § 2.1 for why this is 
so.) The most common contamination is from the host galaxy emission. When a supernova 
is at its brightest this contamination is often small, but in later phases when the supernova 
dims it can become significant and dominant. Host galaxy contamination can be subtracted 
off using spectra of the galaxy when the supernovae is absent, but it is not always clear to 
analyzers of supernova spectra when this has been done. Another kind of contamination 
is peculiar to supernovae associated with gamma-ray bursts (GRBs): such supernovae may 
all be what are now called hypernovae (since they are very energetic and perhaps very 
massive) or as we prefer Type Ic hypernovae since they seem to lack conspicuous lines of 
hydrogen and helium and the strong Si II A6355 line. (Lack of conspicuous hydrogen and 
helium lines together with the lack of a strong Si II A6355 line are the defining observational 
characteristics of Type Ic supernovae.) The contamination is the UVOIR (ultraviolet-optical- 
infrared) afterglow of the GRB itself which tends to be a line-free con tinuum. Correction fo r 



this kind of contamination can be done, but with some uncertainty (IMatheson et al.ll2003l ). 
Yet another contamination (which unfortunately also affects the line pattern) is light echoing 
caused by supernova light from earlier phases reflected off dust clouds and added to the light 
of the current supernova phase. Light echoes, if recognized, can be corrected for, but with 
some error, of course. For the developments in this paper, we assume that contamination of 
all kinds can be adequately corrected for. 

Because of all the effects mentioned above, the continuum shape of supernova spectra 
may be incorrect and misleading about the degree of similarity between spectrum pairs. On 
the other hand the line pattern in the spectra should be better for determining similarities 
since this is less affected by broad-band uncertainties in reddening and calibration. For 
example, if a spectrum pair had an identical line pattern, but different continuum shapes, 
one would have good reason to believe that the supernovae are, in fact, highly similar and 
at a similar phase in their evolution and that the discrepancy in continuum shape was just 
caused by reddening correction and/or calibration error. In the contrasting case, where the 
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line patterns are very different, but the continuum shapes are identical, one would conclude 
the supernovae are not very similar or at different phases in their evolution despite the similar 
continua. In both extreme cases, it is the line patterns that provide the decisive evidence, 
not the continua. Thus, if one can adequately correct or neglect contamination, analysis 
with locally-normalized spectra should give good insight into the intrinsic similarities and 
differences among supernovae. It is a basic premise of this paper that line pattern is a much 
better signature of intrinsic supernova behavior than continuum shape, and so eliminating 
continuum shape information is not too important for spectrum comparisons. 

A second reason for using local normalization applies even when the continuum shapes 
are assumed to be well known. If one is searching just for similarity of line patterns as a 
clue to physical similarities, then varying continuum shapes tend to obscure the similarities 
of the patterns both to the eye and formulaic tests. Being able to measure similarity of 
line patterns for heterogeneous samples of supernovae from varying phases is important in 
studying the time-varying supernova structure. Especially for core-collapse supernovae in 
which even members of the same type show considerably individuality and in which the time 
coverage of their evolution can often be very incomplete, local normalization could become 
an important tool in analysis. 

A third reason, for local normalization is that spectra are frequently analyzed using 
synthetic spectra calculated from highly si mplified radiative transf er: e.g., analyses done 



using the parameterized code SYNOW (e.g-. lBranch et al.ll2003l . l2005l . and references therein) 



Such simplified radiative transfer does not treat continuum radiative transfer with high 
physical realism and in fits of synthetic to observed spectra mismatches in the continuum 
in some regions are obvious and are not considered very important in the important results 
derived from the synthetic spectrum analysis: the important results being line identification 
and ejecta velocity structure. Local normalization of both observed and syn thetic spectra 



can b e used to reduce the distraction of mismatching continuum shapes (e.g., iParrent et al. 



20071 ). 



In § 2 of this paper, we discuss the theoretical basis for local normalization and how 
we actually carry out local normalization. Section 3 gives examples of locally- normalized 
spectra for the Type lib supernova SN 1993J with some analysis of these spectra. The 
UV parts of two of the SN 1993J spectra are Hubble Space Telescope (HST) spectra that 
have not been published before. Also in § 3 (§ 3.5) is an investigation of the continuum 
independence of local normalization. The DIFF1 and DIFF2 formulae are presented in § 4, 
where we also present an example of fitted locally-normalized spectra. As an example of the 
utility of DIFF1/2, § 5 gives some preliminary statistical results for the spectra of supernovae 
of Types lib, lb, Ic, and Ic hypernovae. We will collectively refer to these supernovae as 
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hydrogen-deficient core-collapse (HDCC) supernovae. Conclusions and discussion are given 
in § 6. 

This paper makes use of and refers to material found at the first author's online super- 
nova spectrum database SUSPEND (SUpernovae Spectra PENDing further analysis). The 
URL for SUSPEND is given in a footnote to the abstract. 



2. LOCAL NORMALIZATION 

The basic procedure of local normalization is to divide a spectrum by a smoothed 
version of itself where the smoothing length is sufficiently large that line features are largely 
erased, but the continuum shape is largely unaffected, in the smoothed version. We call the 
smoothing for local-normalization large-scale smoothing to distinguish it from the small-scale 
smoothing used to suppress noise. The spectrum resulting from the division has ideally a 
continuum level of 1 everywhere as judged by the eye. 

In § 2.1, we discuss the theoretical basis for local normalization. In § 2.2, we discuss 
our actual procedure for carrying out local normalization. 



2.1. The Theoretical Basis for Local Normalization 

For concreteness, in our discussion let us model an observed spectrum /A, bs by the 
following heuristic formula: 

/A,obs = S\f\ tCon fx,lm + C\ , (1) 

where S\ is some kind of multiplicative scaling error which could be the effect of unknown 
reddening or calibration error, f\ tCon is the intrinsic continuum spectrum, f\^ n is the line 
spectrum, C\ is some extraneous and unknown contamination flux, and 

f\,int f\,confx,lm j (2) 

is the intrinsic spectrum. The line spectrum is defined by saying it has a height of 1 when it 
is smoothed over some specified smoothing length that is of order the size of the full width 
of an individual line profile: thus 

(/A,lin> = 1 , (3) 

where the angle brackets indicate smoothing with the required smoothing length. We make 
the assumption that S\, f\ tCO n, and C\ are all sufficiently slowly varying with wavelength 
that smoothing them over the smoothing length leaves them effectively unchanged. (We are 
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neglecting the possibility of uncorrected-for light-echo contamination (see § 1).) Now say we 
smooth the observed spectrum over the smoothing length. We obtain 

(fx,obs) = \Sxfx,coaf\,\in) + ^ A ' (4) 

where we have used the assumption that C\ is unaffected by smoothing. Assuming that 
S\j A,con can be approximated as constant over line widths, we now obtain 

{S\f\ t confx,lm) = WA/a,coh) (/ A,lin) = S\f\,con , (5) 

where we have used the assumptions that S\ and /a )C oh are unaffected by smoothing and 
equation (|2D. Making use of equations dlj), dlj), and and our concept of local normaliza- 
tion, the locally-normalized version of /a,ods (which we denote by /a,ln) is given by 

1 + C\/ [S\f\,conf\,lin 



r fx,obs _ S\f\ con f\i[ n + C\ ^ 

JX,LN ~T~f \" c ~f i / A,lin 

\JA,obs/ <JA/A,con + ^A 



l + C x /(S x f x , con ) 



(6) 



What we would really like to have is the line spectrum f\,im ( as we have defined it 
above) and what we can obtain from measurements is /a,ln- The line spectrum f\ t u n is 
what we believe contains much of the information about the object that we are interested 
in. The fx,im and /a,ln spectra are equal if contamination is zero no matter what S\ and 
/A )C on ma y be given our assumptions. (The equality is only exact in the spectrum model of 
this discussion. In reality, it can only be approximate since line and continuum behavior 
cannot be completely separated.) Thus, the effect of any scaling errors are eliminated from 
the locally-normalized spectra in the absence of contamination. This is one of the great 
benefits of using locally-normalized spectra in spectrum analysis. 



If C\/(S\fx t con) « 1 an d /A,iin not much smaller than 1, then 



/a,ln ~ /. 



A.lin 



.con 



(7) 



where we have expanded equation (jSJ) to first order in Cx/(S\f\ iCoa ). We see that contami- 
nation error only cancels to zeroth order. Nevertheless, if C\f (S\fx tCon ) is sufficiently small, 
then contamination will have only a small effect on the locally-normalized spectrum. On the 
other hand, if C\/(S\f\ tCon ) becomes large compared to 1 and /A,ii n , the locally-normalized 
spectrum tends to have value 1 everywhere and be almost independent of the line spectrum 
as one can see from equation (jHJ). 

In the rest of this paper, we assume contamination is, in fact, negligible or has been 
adequately corrected for, and thus that C\ is approximately zero. Given this assumption, 
the locally-normalized spectra are approximately the line spectra as we have defined them 
above. 



- 8- 



2.2. The Procedure of Carrying Out Local Normalization 



There is probably no absolute optimum way of carrying out the basic procedure of 
local normalization. We have developed a prescription, however, that seems very adequate. 
The large-scale smoothing is box-car-like where we march through the spectrum wavelength 
by wavelength and make use of a box-car interval about each wavelength to construct a 
smoothed flux value at that wavelength. The size of the box-car interval is the smoothing 
length. 

The choice of the smoothing length is a key point. If one makes the smoothing length 
too small, line features tend to get suppressed in the locally-normalized spectra: in the limit 
that the smoothing length goes to zero, the locally-normalized spectra become 1 everywhere. 
On the other hand, if one makes the smoothing length too large, then the locally-normalized 
spectra retain continuum features: in the limit that the smoothing length goes to infinity, 
the locally- normalized spectra are just ordinary (or globally) normalized spectra. Unfortu- 
nately in supernova spectra, the wavelength scale of the P Cygni and emission line profiles 
can be very large because of the large Doppler shift velocities in the supernova ejecta. (Su- 
pernova lines in the photospheric or optically-thick phase are broad P-Cygni lines with a 
blueshifted absorption and an emission feature centered roughly on the rest-frame line-center 
wavelength. See, e.g.. IJeffery fc Branch! (11 990 . p. 173-194) for a discussion of P-Cygni line 
formation in supernovae.) Usually, the largest velocities one needs to consider are of or- 
der 30, 000 km s -1 and these lead to relative Doppler shifts of line opacity wavelength of 
v/c ~ 0.1 or of order 10%. However, line features can form at velocities of up to perhaps 



50, 000 km s _1 in Type Ic hypernovae (e.g.. lMazzali et al.ll2000i ) at early times and of up to at 
least 40, OOOkms -1 in other kinds of supernovae if seen at sufficiently early tim es in the UV 
as ev idenced by the P Cygni line of the resonance multiplet Mg II A2797.9 (e.g.. IWiese et al. 



19691 . p. 30; Moore & Merrill 



1968 



, p. 10) seen in the earliest WE spectrum of Type II su- 



pernova SN 1987A (e.g.. IPun et alJll995h . With Doppler shift velocity of 50,000 km 



relative wavelength shift would be v/c ~ 1/6 or of order 17%. Since one has to consider 
opacity shifted both to the blue and the red, the relative wavelength width of a line profile 
in extreme cases could be of order 34 %. Continuum features that vary on the scale of 34 % 
of wavelength certainly exist. Thus, for supernovae at least in extreme cases, one cannot 
completely decouple line and continuum behavior. 

One must make a choice for the smoothing length for locally-normalized spectra that 
does not suppress line behavior too much and that does not let too much continuum behavior 
leak in. There is probably no absolute optimum choice for all cases. We have found that a rel- 
ative smoothing length of 30 % of the current wavelength in the smoothing procedure is good: 
the smoothing region extends 15 % blueward and 15 % redward of the current wavelength 
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as one marches through the wavelength points. The choice of 30 % is motivated as follows . 



The continuum level determined by blackbody-behavior and reddening (ICardelli et al.lll989 



O'Donnelll 119941 ) is relatively constant over a region of 30% surrounding a wavelength. On 
the other hand, as mentioned above, usually the largest velocities one needs to consider are 
of order 30,000kms _1 and these lead to relative Doppler shifts of line opacity wavelength 
of v /c ~ 0.1 or of order 10 %: the whole line profile in this case will vary over ~ 20% of the 
current wavelength. Thus a relative smoothing length of 30 % should be effective at smooth- 
ing away most line features in the smoothed spectrum, but not smoothing the continuum 
shape too much. 

To calculate the large-scale smoothed spectrum itself, a common approach is to simply 
numerically integrate (e.g., using the trapezoid rule) the spectrum over the box-car interval 
centered on the current wavelength and divide by the box-car interval to obtain an aver- 
age flux for the box-car interval which one assigns to the current wavelength. For ordinary 
small-scale smoothing, nothing special needs to be done about the ends of the spectrum 
even though they can be treated somewhat wrongly. As an end is approached, the box-car 
interval gets progressively truncated at that end and the calculated average flux is progres- 
sively a less good approximation for the smoothed flux at the current wavelength: it is the 
smoothed flux for some point in the box-car interval farther from the end. In small-scale 
smoothing, the smoothing length is comparatively short and somewhat bad behavior at the 
spectrum ends is usually unimportant. But with a large smoothing length, such as we re- 
quire for local normalization, this bad behavior at the ends can significantly degrade local 
normalization near the ends. If the spectrum continuum is relatively flat near the ends, there 
may be no noticeable degradation. But if the continuum is rising/falling at an end, then the 
smoothed spectrum becomes too small/large at that end and the locally-normalized spec- 
trum continuum can appear to be rising/falling from 1 at the end. To avoid this problem, 
we fit a line using least-squares to the spectrum in the box-car interval and determine the 
smoothed spectrum at the current wavelength from this fitted line. (Actually we fit a line 
to logarithmic flux as a function of logarithmic wavelength. This means we are fitting the 
flux to a power- law function of wavelength. See below.) In the interior of the spectrum, 
this least-squares line fitting gives almost the same result as calculating an average flux for 
the box-car interval, but at the spectrum ends it usually prevents the locally- normalized 
continuum from departing obviously from 1. Since spectra are not necessarily very linear (or 
like a power-law function) over the box-car interval, it is not immediately obvious that this 
fitting of a line (actually a power-law function in our calculations) should always work well. 
But the resulting locally-normalized spectra seem generally pleasing to the eye including at 
the ends and we accept the eye's judgment for our prescription. 

In plotting spectra, we often prefer to plot logarithmic flux versus logarithmic wave- 



-10- 



length. The choice of logarithmic wavelength owes to the fact that line profile widths are 
determined by large Doppler shifts of line opacity wavelength from the rest-frame line opacity 
wavelength. For example, say a particular component of the line profiles forms in the ejecta 
moving at velocity v along the line-of-sight counting velocity as positive in the direction away 
from the observer. For a given line of rest-frame line-center wavelength Aii ne , the component 
will appear in a spectrum at 



r V I C 

A = Alinc^/^— -j- « A Unc (l + V / c) , (8) 



wher e we have used the relativistic Doppler shift formula (e.g., Lawden 1975 . p. 78; Mihalasl 



19781 . p. 495) and the 1st order (in v/c) Doppler shift formula. (The 1st order Doppler shift 
formula is usually adequate for supernovae.) From equation (JS]), we see for a set of lines that 
the components in a spectrum plot against wavelength will appear at different wavelength 
shifts from the line-center wavelengths. Thus, the line profiles will tend to have varying 
widths that scale with Aii ne . In spectrum figures, this has an unp leasing effect of making 
the blue side of a spectrum look relatively crowded with line features and the red side look 
relatively uncrowded with line features. If one plots using logarithmic wavelength, then the 
aforementioned particular component is found at 

log(A) = log I Ai ine W 7 ' ^J 1 I = l0g( Anne) + log ' ' ' 




v/c I '" ' \V 1 — v/c 

W(l + v/c)} = log(A linc ) + log(l + v/c) « log(A linc ) + (~) log(e) , (9) 

where we have made an expansion of the logarithm to 1st order in v/c which is consistent 
in order of expansion with the 1st order Doppler shift formula. We see that in logarith- 
mic wavelength, the shift of the particular component is the same for all line profiles (from 
relativistic Doppler shift formula) and that this shift is approximately proportional to v/c. 
Thus, the line profiles will tend to have the same width across the spectrum when plotted 
using logarithmic wavelength. Varying line strength and line blending, of course, prevent the 
profiles from having an identical appearance when plotted versus logarithmic wavelength. 
Using logarithmic wavelength gives the plotted spectrum a balanced and pleasing appear- 
ance. The above argument for plotting spectra versus logarithmic wavelength applies for 
both ordinary spectra and locally-normalized spectra. 

In plotting ordinary spectra (i.e., not locally- normalized spectra), we often make the 
choice of plotting logarithmic flux since it gives a more equal appearance to small and large 
features of a spectrum. Say we take equation (2) (see § 2.1) for the intrinsic flux (which for 
the sake of argument we assume that we know) and take the logarithm: we get 

log (/ A ,int) = log (/ A ,con) + log (/ A ,lin) • (10) 
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From equation (10), we see that line profiles of the same relative size (i.e., the same /a,hii 
profile) will have the same absolute size when logarithmic flux is plotted since the line profiles 
are just added to the continuum for logarithmic flux and not multiplied by the continuum 
as in the non-logarithmic case. Of course, varying relative line profile size and line blending 
will make the logarithmic line profiles vary in absolute size. Nevertheless there is some 
equalization in line profile appearance in using logarithmic flux. Why do we prefer the more 
equal appearance of plotting logarithmic flux when large line features, of course, represent 
larger effects? The reason is that small line features present fine tests for spectrum modeling 
and clues about quantities (e.g., composition) that may not be apparent in the large line 
features. Thus, giving more equal weighting to large and small line features makes sense. 

Using local normalization has a similar effect to using logarithmic flux: there is an 
equalization in plotted size of line profiles of the same size relative to the continuum in 
the ordinary spectra. Nevertheless, we plot the logarithmic flux of the locally normalized 
spectra in plotting. This is for consistency with our usual practice with ordinary spectra. 
As mentioned above, in making locally-normalized spectra, we actually use logarithmic flux 
and logarithmic wavelength in creating fitted lines for the large-scale smoothing rather than 
just using flux and wavelength. This procedure in creating the locally-normalized spectra 
tends to make the locally-normalized continuum more like 1 everywhere to the eye on a 
log- log plot. Locally-normalized spectra calculated using a line fitted to flux as a function of 
wavelength (which we will call the unlogged locally-normalized spectra) look similar to those 
calculated with our logarithmic fitting procedure, but the difference is clearly significant to 
the eye. 

Before applying the large-scale smoothing in our local normalization procedure, we also 
apply to all our spectra a small-scale box-car smoothing (using trapezoid rule integration) 
with a relative smoothing length of 5\/\ = 1/300 which corresponds to a Doppler shift 
velocity of about 1000 km s -1 . This smoothing reduces the noise in some of the spectra 
particularly at the spectrum ends. For some particularly noisy spectra, we sometimes use 
a larger small-scale smoothing length. The small-scale smoothing should not degrade the 
spectrum features very much. Supernova line features generally vary relatively slowly over 
wavelength scales corresponding to Doppler shift velocities of order 1000 km s -1 . Interstellar 
and telluric lines can vary rapidly on shorter wavelength scales, but these are not intrinsic 
to supernova spectra and could be eliminated from the spectra if necessary. 

A stand-alone fortran-95 code locnorm.f for making locally normalizing spectra is avail- 
able for downloading from SUSPEND under the heading Useful Programs. All the supernova 
spectra currently in SUSPEND are shown plotted in original and locally-normalized format 
under the heading Supernovae by Epoch in html files in supernova directories. 
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EXAMPLES OF LOCALLY-NORMALIZED SUPERNOVA SPECTRA 



Figures 1, 2, 3, and 4 (which we discuss individually below in separate subsections) show 
locally-normalized supernova spectra plotted on log-log plots together with the original spec- 
tra in the wavelength representation (i.e., f\ or flux-per- unit- wavelength representation): the 
locally-normalized spectra are obviously the ones with continuum levels of about 1 and the 
original spectra have been vertically shifted to be well displayed. The spectra are all for 
Type lib supernova SN 1993 J which occurred in M81 (NGC 3031). T hey have been dered - 
shifted using host galaxy heliocentric velocity — 39±2kms~ 1 from Leda (IPaturel et al.ll2003l ). 
(Since the host galaxy heliocentric velocity is negative, the spectra are, of course, actually 
redshifted to be put in the rest frame.) The original spectra h ave not been dere d dened . The 
E(B — V) reddening value for SN 1993J is quite uncertain. iRichmond et al.l (119941 ) after 
a lengthy consideration of many methods of determining E(B — V) for SN 1993J suggest 
that the SN 1993J E(B — V) value is in the range ~ 0.08-0.32 mag which is a range from 
low to moderate E(B — V). The possible intrinsic spectra obtained by correcting an ob- 
served spectrum using E(B — V) val ues from this range and a standard reddening law (e.g., 
Cardelli et al.lll989l ; IO , Donnelllll994l ) would show a considerable range in continuum behav- 
ior. Thus, the intrinsic continuum behavior of SN 1993J for any phase is quite uncertain. 
Avoiding the uncertainty in dereddening and calibrating observed spectra is one of the main 
reasons for using local normalization as discussed in § [TJ 

Because M81 is a t a distance of only 3.63 ± 0.34 Mpc (as determined from Cepheids 
(IFreedman et al.l Il994l )). SN 1993 J became a very bright and well-obs erved supernova. It 
was discovered on 1993 March 28.86 UT (JD 2449075.36) lRiperolll993n which was probably 
only about a day after expl osion (or core collapse) which may have been 1993 March 27.9 ± 
0.2 (JD 2449074.4 ± 0.2) IRichmond e t al.l |l994h- An initial very hig h peak in UVOIR 
bolometric brightness was partially observed (e.g.. IRichmond et al.lll994l ): this peak, which 
probably happens for all core-collapse supernovae soon after shock break-out, is usually 
unobserved, and so is not used to define conventional UVOIR bolometric maximum light 
even though it may often/always be higher than conventional UVOIR bolometric maximum 
light as it is for SN 1993J. SN 1993J's (conventional) maximum lights occurred on 1993 
April 16 (JD 2449093.7) in B, 1993 April 17 (JD 244 9095.0) in V, and 199 3 April 17 ± 1 
(JD 2449094.7 ± 1.0) in UVOIR bolometric luminosity IRichmond et al.lll994h. The UVOIR 
bolometric luminosity rise time to maximum light was about 20 days (IRichmond et al.lll994l ). 
The maximum B and V apparent magnitudes were, respectively, 11.35±0.05 and 10.80±0.05 
(iRichmond et alill994h . 



SN 1993J was originally classified as a Type II supernova because it had conspicuous 
hydrogen Balmer lines. However, these lines did not become as strong as in typical Type II 
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supernovae and by about 40 days after explos ion SN 1993J evolved to resemble Type lb 



supernovae (e.g., iFilippenko fc Mathesonlll993l ). Type lb supernovae have metal lines and 
conspicuous helium lines and, it seems in general, h ydrogen lines that a re somewhat weak 
and hard to identify especially after maximum light ([Branch et al.ll2002l ). Supernovae that 
undergo a transfor mation from Type II to Type lb have come to be called Type lib's using 
a term invented by IWoosley et al.l (119871 ) in a theoretical prediction that small-hydrogen- 
envelope supernovae could exist. Type lib's are considered to be hydrogen-deficient com- 
pared to ordinary Type II supernovae, but probably have more hydro gen than Type Ib/c's. 



The observational prototype Type lib is SN 1987K (jFilippenkol Il988l ) . Tw o other Type lib 



1999 



supernovae with a sig nificant number of published spectra a re SN 1996cb (jQiu et al 
Matheson et~aDl200lh and SN 1998fa jMatheson et alJl200lh. As of 2007 March 8, the su- 



pernova list of ICentral Bureau for Astronomical Telegrams! (120071 ) contains 25 Type lib or 
possible Type lib supernovae (including SN 1987K which is actually listed as a Type II). 

The ejecta of SN 1993J had a shock interaction with a thick circumst ellar wind shed 



by th e supernova progenitor as evidenced by radio and X-ray emission (e.g., iFransson et al. 



19961 . and references therein). In this respect, SN 1993 J was like th e Type II supernov ae SN 
1979C and SN 1980K which also showed strong radio emission (e.g.. lWeiler et al.lll986l ). The 
interaction also had an effect on the UV spectra of SN 1993J similar to that on SN 1979C 
and SN 1980K. We will discuss this UV effect briefly in § 3.2. 

We have made some line identifications in the figures by labels giving the line's parent 
ion or specific line name in the case of the Balmer lines. The labels are only given for the 
lines in the locally normalized spectra since the corresponding lines in the original spectra 
are identified by comparison to the locally normalized spectra. For clear P Cygni lines, we 
put the identification labels with the blueshifted absorption features of the P Cygni lines 
since this is usually the most distinct part of the line profiles. For emission lines, we put 
the identification labels with the emission features of the emission lines. P Cygni lines 
dominate optical/IR spectra in the photospheric phase and emission lines in the nebular 
phase. Figures 1 and 2 show photospheric phase spectra and Figures 3 and 4 show nebular 
phase spectra. 

We should remark that in the photospheric phase, a supernova is optically thick in the 
optical/IR. As the supernova expands, the opacities must fall with decreasing density and 
the supernova gradually enters the nebular phase where the supernova is optically thin in the 
optical/IR. Supernovae tend to remain optically thick in the UV for much longer than in the 
optical/IR because there are many strong metal lines in the UV particularly from the iron- 
peak elements. There is no sharp dividing time between the two phases. One often just starts 
using the term nebular phase when the emission features start to dominate the optical/IR 
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line pattern. There is also no sharp distinction between P Cygni lines (possessing blueshifted 
absorptions and emission features usually centered on the rest-frame line-center wavelengths) 
and emission lines (which are just emission features usually centered on the rest-frame line- 
center wavelengths). (The centering of emissions on the rest-frame line- wavelengths is often 
approximate and sometimes the emissions are significantly shifted. For example, strongly 
blueshifted emission lines can occur in the UV. For examples of such emission lines, see 
§§ 3.2-3.4.) In practice, when the absorption feature of a P Cygni line has become hard to 
identify, we can say it has become or nearly become an emission line. 

We should also remark that supernova matter after the explosion phase is in uniform 
motion and the whole ejecta structure just scales up linearly with time t since the explosion 
which is effectively a point event. This motion is called homologous expansion. In homolo- 
gous expansion, the radial velocity of any matter element of ejecta is a good comoving-frame 
radial coordinate and we use it as such as is customary in supernova studies. The densities 
of the matter elements in the ejecta scale as t~ 3 . 



3.1. Figure 1: the 1993 March 31 Spectrum 



Because of its closeness and brightness, SN 1993J was spectroscopically very well ob- 
served with the earliest spectrum taken on 1993 March 29.88 UT (JD 24 49076.38) with the 
Isaac Newton Telescope by E. Perez and D. Jones ( Clocchiatti et all 1995h. The spect rum in 
Figure 1 is also an early one from 1993 March 31 (JD 2449078.35) (IBarbon et al.lll995l ) which 
is about 4 days after explosion and is about 16 days before UVOIR bolometric maximum 
light. The original spectrum rises very steeply to the b lue showing that the su pernova pho- 
tosphere was very hot at this phase. From Figure 11 of iRichmond et al.l (119941 ). we estimate 
that a blackbody spectrum fit to the spectrum (if it was well calibrated) for this phase would 
give a blackbody temperature in the range ~ 12000-23000 K where the range is caused by 
the uncertainty in the reddening. The uncertainty in spectrum temperature highlights the 
difficulty in relying on continuum shape for determining the physical properties of super- 
novae from heterogeneous types with uncertain reddening and supports our argument for 
giving more weight to the line pattern than to the continuum shape (see § 1). 



We hay e iden tified the lines in Figure 1 mostly following IBarbon et al.l (119951 ) and 
Baron et al.l (119951 ) . For this spectrum and the others below, we also rely on our gen- 
eral understanding of spectrum formation in making identifications. The B.5 identification 
is uncertain because the absorption feature identifi ed as possibly due to H<5 is rather weak 



and indefinite. The He I line is He I A5875.7 (e.g., IWiese et al.lll966l . p. 14) and the Ca II 



line actually arises from the Ca II IR triplet (i.e., Ca II A8579.1, where we have cited the 
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(//-weighted me an line wavelength : the component lines are at 8498.02 A, 8542.09 A, and 
8662.14A) (e.g., 



Wiese et al.lll969l . p. 251). (Hereafter we usually cite only the (//-weighted 
mean line wavelength for multiplets.) The line velocities (i.e., the velocities correspond- 
ing to the Doppler blueshifts of the line absorption minima from the rest-frame line-center 
wavelengths) are about 11000 km s -1 (H5 if this is the correct identification), 5900 km s™ 1 
(H 7 ), 13400 km s" 1 (H/3), 12800 km s" 1 (He I), 12700 km s" 1 (Ha), and 26800 km s" 1 (Ca II 
IR triplet). We use the ^/-weighted mean wavelengths for calculating the line velocities of 
multiplets. (Here and throughout our discussion, wavelength measurements are made from 
the locally- normalized spectra and we do a little smoothing by eye to determine the minima 
when needed.) The low velocity of H7 may be because of line blending with B.5 since the 
lines are separated in Doppler shift velocity by only 17000 km s _1 . The much higher velocity 
for the Ca II IR triplet line than for the other lines makes the identification a bit uncertain, 
but there seems nothing else it could be. In fact, high-velocity line features, for lines in- 
cluding the Ca II lines, have been found for at least one other HDCC supernova (Type Ibc 



SN 2005bf (Fo lateUi et all 2006 IParrent et all 120071 ): see § □ for HDCC supernovae) and 



the work of IParrent et al.l (120071 ) imples that high- velocity line features might be common 
in Type Ic supernova spectra. (By high-velocity features, we mean line absorptions form- 
ing well above the photosphere at velocities of perhaps > 15000 km s _1 .) Given the above 
remarks, it is possible, speaking speculatively, that HDCC supernovae of all kinds (such as 
Type lib SN 1993J) might have some high-velocity line features in some cases. 

A remarkable feature of the line pattern is that the H/3 line is much stronger than the 
Ha line: this cannot happen in local-thermodynamic equilibrium (LTE). The feature may 
require a non-LTE (NLTE) explanation. On the other hand, the feature may be a result of 
an observational error of some kind since another SN 1993 J spectrum from the same day 

( 19951 ) doe s not clearly show the H/3 line stronger than the Ha 



reported by iBaron et al. 



line. In any case, IBaron et al.l (119951 ) using the NLTE code PHOENIX were able to fit the 
continuum well, but not the lines for the March 30-31 phase of SN 1993J. Their difficulty 
may lie in the inaccurate physical structure of their model and/or an inaccurate reddening 
correction. 



3.2. Figure 2: the 1993 April 15 Spectrum 



Figure 2 shows the SN 1993 J spectrum from 1993 April 15 (jJeffery et al.lll994l ): this 
was 2 days before UVOIR bolometric maximum light and about 18 days after explosion. 
The spectrum is a combination of an HST UV-blue-optical spectrum and Lick Observatory 
ground-based spectrum. The HST spectrum was obtained as part of the HST Supernova 



-16- 



INtensive Study (SINS) General Observer program fcirshner et alJll988h . The two spectra 
agreed well in shape and scale in their overlap region 3120-3276 A, and so we simply joined 
them at 3240 A cut ting o ff the overlapping e nds. The line identifications mostly follow from 



Jeffery et al.l ( 



1994) and iBaron et al.l f 19951 ) . The line velocities are in the range ~ 7000- 



11000 km s 1 (j Jeffery et al.lll994l ). We emphasize though that the lines are often blended and 



the identifications in many cases only indicate the strongest contributor to the line feature. 
The 4 strongest hydrogen Balmer lines, the He I A5875.7 line and the Ca II IR triplet line are 
again present. We note that Ha is now the strongest line in the optical part of the spectrum 
(but note that the Mg II line in the UV is comparably strong: see below) and that it is in 
net emission which is typical of Ha lines in Type II supernovae. 

The absorption centered at about 3835 A i s caused by the Ca II H&K lines: i.e., the 
multiplet Ca II A3945.2 (e.g.. IWiese et al.l Il969l . p. 252). This multiplet is very strong be- 
cause it is a resonance multiplet (i.e., it arises from the ground level). The line velocity for 
the Ca II H&K lines is 8360 km s -1 . The F e II lines mainly aris e from the multiplets Fe II 



A4555, Fe II A4561, and Fe II A5060 (e.g., iKirshner et al.lll993 . Table 4). The absorption 



centered at about 7570 A may be th e absorption of a P Cygni line caused by the multiplet 
O I A7773.4 (e.g., IWiese et al.lll966l p. 79). This line is common in Type Ib /c supernovae 



Branch et al 



20021 ) and Type la supernovae (e.g., IKirshner etaDll993h . but unfortu- 



(e.g., 

nately its absorption usually coincides with a strong telluric line and it is sometimes, as here, 
unclear if the telluric line has been adequately corrected for. We regard this identification 
as tentative. If the identification is correct, the O I A7773.4 line velocity is 7830 km s -1 . By 
about 4 days past UVOIR bolometric maximum light, the O I A7773.4 line seems clearly 
present (IBarbon et al.lll995l ). 



The UV of SN 1993 J is similar to that of SN 1979C and SN 1980K and is probably greatly 
modified from what a bare supernova would have given by c i rcumst ellar sh ock interaction 
( Jeffery et al.lll994l ). In a series of papers, i Franssonl Jl982l . Il984al lbh and iFransson et aJJ 



(Il984d ) have given a detailed model of supernova UV behavior with circumstellar interaction: 
this work was stimulated by International Ultraviolet Explorer (IUE) data for SN 1979C 
and SN 1980K. In brief, some fraction of UVOIR light from the supernova photosphere is 
scattered by shock-heated electrons in the shock region of ejecta and circumstellar matter 
and Comptonized (i.e., blueshifted up to the UV or extreme UV). Also there will be X-ray 
emission from the shock-heated gas. Some of the Comptonized light and X-ray emission 
just escapes the supernova environment, but some back-heats the outer supernova ejecta 
and creates a complex, layered, photoionization region in the outer supernova atmosphere. 
The emission from back-heating adds to the direct emission from the supernova interior that 
emerges at the photosphere. In the optical, this emission from back-heating is probably 
negligible (except in some lines), but in UV blueward of perhaps ~ 2800 A it appears to 
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become imp ortant by compari son to spectra from supernovae with little or no circumstellar 



interaction (jJeffery et al.lll994l ). As well as continuum emission from recombination and free- 
free emission, UV emission lines from ionized atoms can make a considerable contribution: 
the prominent UV emission lines are primarily resonance lines that have been collisionally 
excited by the back-heated medium. The layers of Comptonized-flux photonization form a 
geometrically narrow shell of ejecta moving at velocities of order 10 kms -1 . The observer 
sees the emission from back-heating primarily from this shell. Below the shell is the optically 
thick photosphere. Thus, the observer mostly sees the emission from back-heating from the 
near-side of the ejecta with line-of-sight velocities ranging from much less than 10 4 kms _1 
away from the observer up to line-of-sight velocities of order 10 4 kms _1 toward the observer. 
For UV emission lines, the varying Doppler shift given by the moving shell will result in 
broad blueshifted emission features. 

An effect that tends to enhance blueshifted line emission from a shell surrounding an 
opaque photo sphere is diff use reflection toward the observer of the line emission from the 



photosphere (jChugail ll988). This reflected line emission must mainly come from the near 
side of the photosphere, and thus will tend to be mainly blueshifted from the rest-frame 
line-center wavelength and hence the enhancement of the blueshifted line emission. If the 
interior of the shell were transparent, the light that could otherwise be reflected would simply 
stream away from the observer and all redshifted line emission from the far side of the shell 
would come toward the observer. The reflection can be caused by scattering from electrons, 
but in the UV, particularly in the nebular phase when the electron optical depth through 
the ejecta has become small, the reflection is probably mainly by scattering from the many 
iron-peak-element lines in the UV. To understand the reflection effect clearly, imagine an 
ideal case of a perfectly reflecting photosphere just interior to a very geometrically thin 
line-emitting shell. Virtually all the line emission from the near side of the shell moving 
inward would be reflected from the near side of the photosphere and some would go to the 
observer. Virtually no emission or reflection would come to the observer from, respectively, 
the far side of the shell or the far side of the photosphere. Since the photosphere matter is 
expanding homologously, the reflected light toward the observer will be mostly blueshifted. 
Note that in moving inward toward the photosphere, the line emitted photons will redshift 
from the rest-frame line-center wavelength in the comoving frame and that the photosphere 
is moving a little more slowly than the shell. (Propagating photons always redshift in the 
comoving frame in homologously expanding atmospheres no matter what the degree of special 



relativistic effects (jJefferyill993l ).) Consequently, the light reflected toward the observer will 
tend on average to be less blueshifted in the observer frame than the line emission coming 
directly from the line-emitting shell. Because the photosphere is just interior to the shell in 
our ideal case, this redshift effect will be small. In a real case, the shell will not necessarily be 
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very geometrically thin, and so some significant amount of redshifted flux from the far side of 
the line-emitting shell could reach the observer. Also some fraction of the observer-directed 
reflected light (mainly from near the photosphere limb where the line-of-sight velocity of the 
photosphere is small) will end up a bit redshifted from the line-center wavelength in the 
observer frame because of the redshift in the comoving frame that happens while the light 
moves inward toward the photosphere. But the contribution of this redshifted reflected light 
is probably relatively small in most cases. Note also, a real photosphere will not, of course, 
be perfectly reflecting: some of the line emission will be absorbed and its energy dispersed 
to other parts of the spectrum. 

We have not attempted a full identification of features seen in the UV region of the 
spectra we present. However, the strong emission-like feature peaking at 2717 A in the 



spectrum must be the resonance multiplet Mg II A2797.9 (e.g., IWiese et al.l Il969l . p. 30; 



Moore fc Merrilll Il968l . p. 10) blueshifted and in emission and produced as discussed by 
Fransson. (The emission features most fully analyzed by Fransson were primarily those from 
blueward of 2000 A, but his a nalysis seems to apply to Mg II A2797.9 as well.) In SN 1979C 
UV IUE spectra (1lNESll2000l ). a similar blueshifted Mg II A27 97.9 emission appea red possi- 
bly about 7 days after optical maximum light (1979 April 22) JPanagia et alill98oh . became 
definite by 33 days past optical maximum (1979 May 18), and was still clearly present in 
the last IUE observation about 111 days past maxim um (1979 August 4 ). (Optical max- 
imum light for S N 1979C wa s about 1979 April 15 (fearbon et alill982h .) In SN 1980K 
UV IUE spectra (jlNESll2000l ). a similar blueshifted Mg II A2797.9 emission seems to have 
arisen about 26 days after optical maximum light (1980 October 30) and be present on the 
last /[/i? observation from about 62 days past optical maximum (1981 January 5). In the 
SN 1980K spectra, the emission is generally less clear than in SN 1979 spectra, but it defi- 
nitely present about 35 days past optical maxi mum (1980 December 9). (Optical maximum 
light of SN 1980K was about 1980 November 4 (|Butalll982h .) We note that an Mg II A2797.9 
emission line is a common feature of quasar emission-line clouds (e.g., iRei chard et al.l 120031 ) 
which may somew hat resemble the UV emission shell of circumstellar-interacting supernovae 
(jFranssonl 1 1 984bl ) . 



The Mg II A2797.9 emission peak wavelength 2717 A corresponds to a Doppler blueshift 
velocity of 8700 km s - 1 which is within the range of P-Cygni line velocities in the optical 
~ 7000-11000 km s" 1 flJefferv et al.lll994h . If we just define 2700 A as the characteristic blue 
edge of the Mg II A2797.9 emission (which is good for later phases: see §§ 3.3 and 3.4 below), 
the corresponding Doppler blueshift velocity is 10500 km s" 1 . We take this velocity as the 
characteristic velocity of the Mg II A2797.9 emission shell. 



Without strong circumstellar interaction, the Mg II A2797.9 multiplet would probably 
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contribute along with Fe II lines to a blueshifted P-Cygni absorption feature. For Type la 
supern ovae, the absorption can be seen in the HST spectra of SN 1992A in the photospheric 



phase (IKirshner et al. 



19931 ) and, perhaps, well into the nebular phas e at 291 days after B 



maximum light although not so identified ( Ruiz-Lapuente et al. 19951). For Type II super- 
novae, the absorption can be seen in the IUE spectra of SN 1987A flPun et al.lll995h . The 
difference for the Mg II A2797.9 multiplet between having blueshifted emission and having 
blueshifted absorption is probably explained by saying that with circumstellar interaction 
the line-forming region is hotter than the inner ejecta and without circumstellar interaction 
it is colder. 

No other identifications in the UV of SN 1993J are probably possible without a more 
detailed analysis. We leave that to future work sine die. 



3.3. Figure 3: the 1993 September 15.5 Spectrum 



Figure 3 shows the SN 1993J spectrum that we have constructed from a 1993 Septem- 
ber 17 HST UV-blue-optical spectrum and a 1993 September 14 ground-based spectrum 
taken at the MMT Observatory by P. Challis and C. S. J. Pun. We will call the com- 
bined spectrum the 1993 September 15.5 spectrum: thus it comes from about 151 days after 
UVOIR bolometric maximum light and 171 days after explosion. Both component spectra 
were obtained as part of the SINS program and neither have been published before (j Challis 
20061 . for the HST spectrum; iBlondinl 120071 . for the MMT sp ectrum). The HST spectrum 
has not been given a definitive flux calibration (jChalliall994l ). The two spectra agreed well 
in shape in the region 4250-4781 A, and so we simply and joined them at 4253.6 A cutting 
off the overlapping ends. The ground-based spectrum had to be scaled by a factor of 0.35. 

The supernova at the time of this spectrum is clearly in the nebular phase because 
of the strong optical/near-IR emission lines. The line identific ations for the 3 st r ongest 
emission lines in the optical/near-IR are based on analysis of iHouck fe Franssonl (119961 ) 
for a 1993 August 15 SN 1993J spectrum: the extrapolation of the identifications to 1993 
September 15.5 is reliable because these lines are strong, persistent, and well identified in 
other HDCC supernovae. The actua l lines are the for bidden multiplet [O I] AA6300, 6364 
which is a resonance mult iplet fe.g.. Iwiese et aL 1966 . p. 82) and the forbidden multiplet 



[Ca II] AA7291,7324 (e.g.. IWiese et al 



1969 



p. 255) and, again, the Ca II IR triplet. 



The absorption with minimum at 3829 A is caused by Ca II H&K lines (the resonance 
multiplet Ca II A3945.2). The line velocity of the Ca II H&K lines is 8800 km s" 1 . There is 
no obvious strong emission feature associated with these lines. Evidently, the Ca II H&K 
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lines are optically thick and strongly trap photons in the spatial resonance regions of the 
lines. The trapping process is as follows. As photons are propagating through the supernova 
atmosphere, they redshift into resonance with the lines in some location (i.e., a resonance 
region) and are absorbed and then resonantly emitted, but cannot easily redshift out of 
resonance because of continually being re-absorbed and re-emitted in the lines after traveling 
and redshifting very little. So there are many re-absorptions and re-emissions. Much of the 
photon energy ultimately escapes by being emitted b y the Ca II IR triplet which shares 
the same upper level with the Ca II H&K lines (e.g., iMoore fc Merrill! Il968l . p. 12): this 
energy leak mostly suppresses the emission feature of the Ca II H&K lines. The lower level 
of the Ca II IR tri plet is metastable and the upper level of the forbidden multiplet [Ca II] 



AA7291,7324 (e.g.. lMoore fc Merrill! Il968l . p. 12). Collisions are weak in the nebular phase 
due to low density, and hence the transitions from the metastable lower level of the Ca II 
IR triplet to the ground level become increasingly radiative which results in the strong [Ca 
II] AA7291, 7324 emission feature. The process of energy from absorption in t he Ca II HfcK 



lines being emitted in other redder lines is an example of fluorescence (e.g., lMihalaslll978l . 
p. 22). 

In the UV, the blueshifted Mg II A2797.9 emission line now peaks at 2726 A which 
corresponds to Doppler blueshift velocity of 7700kms _1 and has now increased in height to 
about 3.4 times the continuum level. Its overall shape, however, is much the same as earlier 
on 1993 April 15 (see Fig. 2 in § 3.2 and Fig. 5 in § 3.4). 

The peak centered at 4563 A may be at least in p art th e emission feature of the semi- 



forbi dden line Mg I] A4571 fe.g.. lMoore fc Mer rill 



196 8i. p. 16: IWiese et al.lll969l . p. 26: lNIST 



20071 ). The emission is relatively weak, but lHouck fc Franssonl (119961 ) predict Mg I] A4571 as a 



weak emission line (seemingly blended with some other emission line) for 1993 August 15 and 
believe the Mg I] A4571 emission should become stronger with time. If the 4563 A emission is 
primarily Mg I] A4571, then the line has been Doppler blueshifted by velocity ~ 500 km s -1 . 
The 1993 November 15.5 spectrum in Figure 4 (see § 3.4) shows a stronger emission peaking 
at 4551 A and if this is caused by Mg I] A4571, then it has been Doppler blueshifted by 
velocity ~ 1300 km s -1 . The 1993 November 15.5 spectrum emission peaking at 4551 A, 
however, has a secondary peak at almost exactly 4571 A. This suggests that the main peak 
is caused by some other line in emission and not Mg I] A4571. In order for the 4563 A and 
4551 A emissions to be primarily Mg I] A4571, there must be some blueshifting asymmetry 
in the ejecta or some significant continuum optical depth. (A significant continuum optical 
depth would tend to suppress the red side of an emission line.) We note that the [O I] A6300 
emission peaks at 6302 A (corresponding to a Doppler redshift velocity of ~ 100 km s -1 ) in 
the 1993 September 15.5 spectrum and at 6300 A (corresponding to virtually no Doppler shift 
at all) in the 1993 November 15.5 spectrum. The lack of significant shifts in the [O I] A6300 
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emission peaks argues for symmetric ejecta and optically thin ejecta in the optical. There 
is the complication, however, that the blue side of the optical generally has more significant 
weak lines than the red side of the optical and may persist being somewhat optically thick 
longer than the red side due to a quasi-continuum opacity of the weak lines. The upshot of 
all these considerations is that we draw no conclusion about the identities of the 4563 A and 
4551 A emissions. These emissions may be blueshifted Mg I] A4571 or Mg I] A4571 blended 
with some other emission line, but more analysis is needed. 

(The case for Mg I] A4571 identification in the optical is different than for the Mg II 
A2797.9 identification in the UV. We expect the nebular-phase ejecta to be optically thin or 
relatively optically thin in the optical and optically thick in the UV at least for a long time 
into the nebular phase. The UV optical thickness in the nebular phase is the cause of the 
blueshift of the Mg II A2797.9 emission as we argued in § 3.2.) 

We note that the strong emission features in the optical/near-IR in both versions of 
the spectrum, but particularly in the locally-normalized version, seem to be adjacent to 
absorptions. There may be real strong absorptions at these locations from lines that are 
optically thick, but a synthetic spectrum analysis is probably needed to show that. The 
locally-normalized spectra for the nebular phase, however, may be a bit deceptive. The 
process of local normalization will tend to make emission features sit in troughs that are 
below 1, and thus create the appearance of adjacent absorptions. In reality, the emission 
lines may just be added to a fairly smooth continuum that holds across the spectrum (or 
at least the optical/IR spectrum) and there may be no significant absorptions (at least in 
the optical/IR spectrum). So in regard to nebular spectra, local normalization as we have 
implemented it may not be ideal in that it probably cannot scale a real physical nebular 
continuum to nearly 1 at least in the vicinity of strong emission features. On the other 
hand, there may be no well-defined physical continuum level in which case applying local 
normalization does not obviously worsen the appearance of the spectra for analysis. Our 
arguments for local normalization based on avoiding dealing with uncertain reddening and 
calibration error (§ [I]) are still valid for the nebular phase. Figure 3 and Figure 4 (§ 3.4) do 
show that local normalization even in the nebular phase will, of course, tend to cause the 
average flux level to be near 1. 

For the photospheric phase, local normalization is probably much better at setting the 
phyisical continuum to near 1 than for the nebular phase. First, there is a significant real 
physical continuum in the spectrum (although it is probably almost never a pure blackbody 
continuum) which is not necessarily the case in the nebular phase. Second, the smoothing 
of the P Cygni absorptions and emissions to create the large-scale smoothed spectrum is an 
averaging that must put the large-scale smoothed spectrum close to the physical continuum 



-22 - 



level. The locally- normalized spectrum will then have a continuum close to 1. 



3.4. Figure 4: the 1993 November 15.5 Spectrum 



Figure 4 shows the SN 1993J spectrum that we have constructed from a 1993 Novem- 
ber 14 HST UV-blue-optical spectrum and a 1993 November 17 Lick Observatory ground- 
based spectrum both of which w ere obta i ned a s part of the SINS program. The HST spectrum 
has not been published b efore (jChallisI 120061) ; the Lick spectrum was taken by A.V. Filip- 
penko and T. Matheson f Filippenkol 1994 ). We will call the combined spectrum the 1993 
November 15.5 spectrum: thus it comes from about 212 days after UVOIR bolometric max- 
imum light and 232 days after explosion. The HST spectrum has not been given a definitive 
flux calibration (jChallislll994l ) and the groun d-based spectrum is not absolutely calibrated, 
although its relative calibration is excellent (jFilippenkol 1 19941 ) . The two spectra agree well 
in shape in the overlap region 3120-4781 A, and so we simply joined them at 4200 A cutting 
off the overlapping ends. The ground-based spectrum, after converting it from the frequency 
representation to the wavelength representation, had to be scaled by a factor of 5.55 x 10~ 9 
to account for overall calibration differences and for unit transformation. The HST spectrum 
blueward of ~ 2425 A seems to decline in a manner too steep to be physically real and we 
do not trust it (nor the locally-normalized spectrum we derive from it) there. 

Qualitatively, the 1993 November 15.5 spectrum is much like the 1993 September 15.5 
spectrum. (We are relying on the locally-normalized spectra for comparisons.) The Mg II 
A2797.9, [O I] AA6300, 6364, and [Ca II] AA7291, 7324 emissions have increased in height 
relative to the continuum and the Ca II IR triplet emission has decreased. The peak of 
the Mg II A2797.9 emission is at 2726 A just as in the 1993 September 15.5 spectrum. The 
emission at peaking at 4563 A (whatever it may be) in the 1993 September 15.5 spectrum 
has increased and shifted so that it now peaks at 4551 A. (See the discussion of this emission 
feature in § 3.3.) The Ca II H&K absorption is not as deep as in the 1993 September 15.5 
spectrum. 

In Figure 5, we show the time development of Mg II A2797.9 emission in locally- 
normalized spectra. Although the emission line scales up with time, its appearance is 
remarkably stable. Our definition of 2700 A as the characteristic blue edge of the Mg II 
A2797.9 emission (with corresponding Doppler blueshift velocity 10500 km s _1 ; see § 3.2) is 
valid for all the phases we display. The behavior of the spectra in the region ~ 2780-2820 A 
is not very certain since we have just omitted wavelength points that correspond to the 
interstellar Mg II A2797.9 lines which arise from gas clouds along the line of sight to the 
supernova. Such clouds, identified from interstellar lines of various ions, have heliocentric 
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velo c ities in the range from about —140 to 181 km s 1 (jWamsteker et al.lll993t IWheeler et al. 
1993h . 



Later HSTXJV spectra of SN 1993 J from day +649 to about day +2561 relative UVOIR 
bolometric maximum light show the Mg II A2797.9 emission as box-like in appearance and 
more symmetrical abo ut the rest-frame line-center wavelength than in the spectra we show 
( iFransson et al.l 120051 ) . The Mg II A2797.9 emission in the later spectra have a Doppler 
velocity full-width of ~ 17000 km s" 1 . It seems that the circumstellar interaction back- 
heating persists through the later phases, but the optical depth through the ejecta has 
greatly diminished, and thus occultation of the far side of the Mg II A2797.9 shell has 
greatly decreased. The flattish top of the late Mg II A2797.9 emission line is charac teristic 
of emission lines from expanding, optically thin shells as first s hown by iBealsH 19311 ): other 
references for the flat-t op emission from expanding shells are iMihalad (119781 . p. 477) and 
Jefferv & Branch! <jl99oL p. 190). 



3.5. The Continuum Independence of Local Normalization 

The main goal of local normalization is to eliminate variations in the continuum in 
comparing spectra. Implicit in this goal is the requirement that locally-normalized spectra 
be virtually independent of the continuum shape. For example, if this were not the case, 
intrinsic spectra that were identical, but with observed continuums distorted in various ways, 
could be transformed to different locally-normalized spectra and one would reach the false 
conclusion that the intrinsic spectra were not identical. To test local normalization for 
continuum independence, we have weighted each of the four original spectra in this section 
by A 2 (effectively converting the spectra to the /„ representation) and A 4 and then locally 
normalized these weighted spectra and compared them to each other and the unweighted 
locally-normalized spectra. 

For each original spectrum, a plot (using an expanded scale compared to that spectrum's 
plot in Fig. 1-4) of the corresponding locally-normalized spectra calculated with and without 
weights shows exact overlap to the eye of the locally-normalized spectra. On the other 
hand, the unlogged locally- normalized spectrum (see § 2.2) corresponding to an original 
spectrum was distinct from the other locally-normalized spectra for that original spectrum 
when plotted: the unlogged locally-normalized spectrum was generally a bit lower than the 
other locally-normalized spectra. 

As mentioned in the Introduction (§ 1), DIFF1 is essentially a mean relative difference 
between the line patterns of two spectra, and thus the reader can roughly understand DIFF1 
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values before seeing the precise formula (which is given below in § 4.1). The DIFF1 values 
for spectrum pairs consisting of the locally-normalized weighted spectra and the locally- 
normalized unweighted spectra (all derived from the same original spectrum) were of order 
0.001-0.0015 for weight A 2 and of order 0.002-0.003 for weight A 4 . These are very small 
mean relative differences. On the other hand, DIFF1 tests of the unlogged locally-normalized 
spectra with the corresponding other locally-normalized spectra gave DIFF1 values of order 
0.4-0.7 which are comparable to DIFF1 values between locally-normalized spectra of the 
same supernova from phases differing by several days. 

We conclude that local normalization will transform continuum-distorted versions of an 
original spectrum to virtually the same locally-normalized spectrum. Thus, local normal- 
ization is sufficiently continuum-independent to lead to valid results in comparing locally 
normalized spectra. On the other hand, different local normalization procedures can pro- 
duce noticeably distinct locally-normalized spectra for a given original spectrum. Thus, one 
should choose a single local normalization procedure and stick with it for valid comparisons 
of locally-normalized spectra. 



4. DIFF1 AND DIFF2 

In this section, we introduce the formulae for DIFF1 (§ 4.1) and DIFF2 (§ 4.2), show 
an example of two locally-normalized spectra fitted by minimizing the DIFF2 function value 
(§ 4.3), and compare DIFF1 and DIFF2 (§ 4.4). 



4.1. DIFF1 

Say fa is the locally-normalized flux at wavelength Aj. Let 

Si = fi-1 (11) 

be the flux difference at wavelength A«. The flux difference is the relative difference of the 
line pattern from the continuum for the original spectrum or the absolute difference of the 
line pattern for the locally-normalized spectrum. We define DIFF1 between a spectrum 1 
and a spectrum 2 by the formula 

DIFFl = iV r , (12) 

1^ max (|5 M |, |^|,0 1 ' 

where the subscripts 1 and 2 indicate the spectrum and / is the total number of wavelength 
points in the summation. (The wavelength points are equally spaced in logarithmic wave- 
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length: see below.) The £ is a small number that we set to 10" 15 . This is typically the 
smallest number that added to 1 that creates a number significantly different from 1 with 
fortran 95 double precision arithmetic. The £ prevents arithmetic exceptions, but is only 
very rarely invoked. 

From equation (12), it is clear why we describe DIFF1 as the mean relative difference 
in line pattern between two spectra (see §§ 1 and 3.5). If DIFF1 << 1, the spectra will be 
much alike. If DIFF1 > 1, the spectra will be substantially different. For example, if the 
flux differences (i.e., the <Vs) for one spectrum in a DIFF1 test were almost all only a small 
fraction of those for the other spectrum (i.e., one spectrum had much weaker lines than the 
other), the DIFF1 value would be close to 1. In fact, the DIFF1 function has an upper limit: 

\S±,i - 5 2 ,i\ < \5i,i\ + \S 2 ,i\ < 2max(|5 M |, |<J 2l *|) < 2max(|5 M |, \6 2 ,i\,Q , (13) 

and thus 

DIFF1 < 2 . (14) 

We would not usually expect a DIFF1 value approaching 2 since that would mean that two 
input locally-normalized spectra had line patterns that were nearly mirror images of each 
other about the continuum level: i.e., the line patterns would be strongly anticorrelated. In 
fact, any DIFF1 values well over 1 would show some anticorrelation between spectrum line 
patterns. There is no strong physical reason to expect much anticorrelation of supernova 
spectrum line patterns, but some accidental anticorrelation should happen, and so some 
DIFF1 values over well over 1 may occur when applying DIFF1 to heterogeneous samples of 
spectra. 

We have subjected the set of all possible spectrum pairs drawn from the sample of HDCC 
supernovae described in § 5 to the DIFF1 test. (There are 17 supernovae, 168 spectra, and 
27588 valid spectrum pair tests: some possible pair tests are excluded as invalid by the 
validity criterion discussed below.) The resulting distribution of DIFF1 values has mean, 
standard deviation, minimum value, and maximum value of, respectively, 0.865, 0.137, 0.244, 
and 1.320 (where we have reported more digits than are significant to allow for numerical 
consistency checks). That the mean 0.865 is so close to 1 (actually within la of 1) shows that 
spectrum pairs on average are not highly similar. This is understandable since the spectra 
come from a heterogeneous sample of supernovae and from many phases: the spectrum of 
a supernova can evolve significantly with phase. The standard deviation value shows that 
spectrum pairs significantly (i.e., la) more alike than average will have DIFF1 values of 
< 0.73. The minimum 0.244 suggests that even spectra from the same supernova at close 
to the same phase vary from each other significantly. Note we have excluded redundant 
spectra from the same phase period (which we usually set to being 1 day) for a supernova 
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(as described in § 5.1) since those should be nearly identical to the included spectrum aside 
from observational error. As we expected, the maximum DIFF1 value is much less than 2. 

We note that if we used /j's in the formula for DIFF1 (i.e., eq. (12)) instead of S^s 
(which would only change denominators, in fact), we would have lessened the sensitivity 
to line patterns. For example, Type lib and Type Ic hypernovae at early phases (i.e., well 
before maximum light) both have relatively weak lines in the optical. In the case of Type lib 
supernovae, this is because they have hydrogen-dominated atmospheres that when hot (i.e., 
over of order 10000 K) show weak lines. (Why the lines should be as weak as they are in 
early Type lib sp ectra has not yet been fully elucidated theoretically by NLTE analysis 



(IBaron et al.lll995l ).) Type Ic hypernovae show weak lines at early times because they are 
hot (which tends to cause weak optical lines also in cases with little or no hydrogen) and 
because their velocity structure is very fast giving them ext remely broad P-Cygni lines with 



line v elocities of up to 30000 km s 1 and perhaps higher (e.g. , iMazzali et al.ll2000i ; iPatat et al. 



200ll ). The stretching out of P-Cygni lines by very high velocities tends to make them 
shallower and more overlapping. Both effects tend to give line features a small vertical scale 
relative to the continuum. If we used f^s in the formula for DIFF1 instead of 5{S, the 
DIFF1 values between early phase spectra of Type lib's and Type Ic hypernovae would be 
small even though the spectra and the supernovae are intrinsically quite different because 
the formula numerators would be relatively small and the formula denominators relatively 
large. However, the given formula for DIFF1 distinguishes the two types of spectra because 
the formula denominators become small when the lines features are small. 

For our calculations, we interpolated the locally-normalized spectra onto a grid of 
equally spaced points in (natural) logarithmic wavelength. For coding convenience, we chose 
ln(l A) = to be the zero point of the logarithmic wavelength scale. For the grid increment, 
we chose 

Aln(A) = . (15) 

v ' 3000 v ; 

If the grid wavelengths labeled by % are Aj and we define AAj = \ — Aj_i, we find that 

A ln(A) = In(Ai) - ln(Aj_i) = In ( -^-) =\n(l + « ^ (16) 



_ Aj-i / \ Aj_i / A 

to first order in small AAj/Aj_i. We now see that the logarithmic wavelength increment 
corresponds to a relative change in wavelength for each increment of approximately 1/3000. 
Thus, near 3000 A, the wavelength increment is about lA; near 6000 A, about 2 A. The 
Doppler shift velocity corresponding to a relative wavelength increment of 1/3000 is about 
lOOkrns -1 . As mentioned in § 2.2, supernova line features generally change relatively slowly 
over wavelength shifts corresponding to velocity changes of order lOOOkms -1 , and thus our 
grid should be sufficiently fine for almost all supernova spectra. As also mentioned in § 2.2, 
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interstellar and telluric lines can vary on shorter wavelength scales, but these are not intrinsic 
to supernova spectra and could be corrected for if necessary. 

We only apply the DIFF1 test to what we call the good overlap region in logarithmic 
wavelength of two locally-normalized spectra. In this logarithmic wavelength region, both 
the spectra seem physically good. By physically good, we mean not excessively noisy and 
not showing unphysical behavior in a narrow wavelength band because of some calibration 
error. We only allow one good region for a spectrum: i.e., we never multiple good regions 
for a spectrum which are separated in wavelength. This is not much of a limitation since the 
bad regions are usually the spectrum ends which can often be very noisy and have narrow 
wavelength band calibration errors. We only consider DIFF1 tests valid if the good overlap 
region in logarithmic wavelength is > 0.2 which corresponds to a good overlap region in 
relative wavelength of ~ 22 % relative to the lower wavelength of the region. Since almost 
all the spectra we consider are good or nearly good over the range 4000-7000 A (logarithmic 
wavelength range ~ 0.5596 or relative wavelength range of 75% relative to 4000 A), almost 
all DIFF1 evaluations are considered valid. 



The only difference from equation (fl2~i) for DIFF1 (§ 4.1) is the term k in the wavelength 
subscript of one of the flux differences (defined by eq. (Hip in § 4.1). This k is chosen to 
minimize the DIFF2 function value: this minimized value is the DIFF2 value itself. We 
allow up and down (natural) logarithmic wavelength shifts of only up to 0.05. Thus, we 
allow up and down relative wavelength shifts of only about 5 %. Given that the logarithmic 
wavelength grid increment is Aln(A) = 1/3000 (see § 4.1, eq. ffTBT) ). it follows that we allow 
k to vary up and down from zero by 0.05/(1/3000) = 150. If the minimizing value of k turns 
out to be one of the k limits —150 or 150, then we consider the test invalid since the true 
minimizing value of k is almost certainly beyond the k limits. This occasionally happens 
and indicates the tested spectra are likely not very similar and likely cannot be made to look 
similar by shifting them relative to each other in logarithmic wavelength: if they can, the 
resemblance is likely accidental. Also we only consider DIFF2 tests valid if the good overlap 
region in logarithmic wavelength (which is specified in § 4.1) when the locally- normalized 
spectra are not relatively shifted (i.e., when k = 0) is > 0.2 which corresponds to a good 
overlap region in relative wavelength of ~ 22 % relative to the lower wavelength of the region. 



4.2. 



DIFF2 



The DIFF2 formula is 




(17) 
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This is the same overlap criterion for validity as used for the DIFF1 test and, as mentioned 
in § 4.1, is almost always met for spectra we consider. 

By its nature, DIFF2 < DIFF1 for any spectrum pair. Thus, DIFF2 has an upper limit 
of 2 like DIFF1 (see § 4.1). But it is even less likely for DIFF2 than for DIFF1 that any 
value will approach 2. To do that would imply that the locally-normalized spectra were 
anticorrelated no matter how one shifted them over the allowed shifting range. Since even 
accidental anticorrelations of spectra in a large sample will tend to be avoided by the shifting, 
we do not expect DIFF2 values that are well over 1. 

We have subjected the set of all possible spectrum pairs drawn from the sample of 
HDCC supernovae described in § 5 to the DIFF2 test. (There are 17 supernovae, 168 spec- 
tra, and 25172 valid spectrum pair tests: some possible pair tests are excluded as invalid by 
the validity criteria discussed above.) The resulting distribution of DIFF2 values has mean, 
standard deviation, minimum value, and maximum value of, respectively, 0.784, 0.116, 0.229, 
and 1.131 (where we have reported more digits than are significant to allow for numerical 
consistency checks). The mean, minimum, and maximum values are, of course, a bit smaller 
than the corresponding values for DIFF1 reported in § 4.1: all the DIFF2 values that go into 
determining these distribution characteristics are equal to or smaller than the corresponding 
DIFF1 values used to evaluate the distribution characteristics in § 4.1. The standard de- 
viation value shows that spectrum pairs significantly (i.e., la) more alike than average will 
have DIFF2 values of < 0.67. The maximum value as expected is not well over 1. 

The rationale for DIFF2 is as follows. The structure of P-Cygni line profiles in su- 
pernovae is correlated with the velocity structure of the ejecta since that velocity structure 
determines the Doppler shifts of line features. Most noticeably, weak, unblended lines in the 
photospheric phase tend to have their absorption featu res reach a minimum at a wavelength 



shift corresponding to the photospheric velocity (e.g.. iJeffery fc Branch! Il990l . p. 188). All 
parts of the absorption feature have blueshifts that are usually correlated with the pho- 
tospheric velocity in the sense that the greater the photospheric velocity, the greater the 
blueshift of each part of the absorption feature: this is just because the whole velocity struc- 
ture of the line-forming region at any phase tends to be correlated with the photospheric 
velocity at that phase. In comparing photospheric spectra for similar supernovae, there is 
often some difference in the width scale of line profiles, particularly the absorption features, 
attributable to some difference in the velocity structure. The difference may be intrinsic to 
the supernovae: one supernova may just have layers that move faster than the corresponding 
layers in the other supernova at the same phase. On the other hand, the difference may just 
be a matter of phase since the photosphere recedes into the ejecta as time passes because 
the overall optical depth of the supernova falls. Of course, intrinsic and phase differences in 
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velocity structure are usually both present to one degree or another. 

The eye can usually discern similar line profile patterns despite differing velocity struc- 
ture. However, the DIFF1 test may not pick out a striking similarity in two spectra seen by 
the eye because Doppler-shift-caused offsets of line profiles between the two spectra may lead 
to large contributions to the DIFF1 value. One can imagine compensating for the different 
velocity structures by a wavelength varying shift, but this seems too complex for practical 
use. Since the blueshifted absorption features are often larger than the emissi on features 



(necessarily the case if the lines are pure scattering (e.g.. iJeffery fc Branch! 1 199(J . p. 189)), a 
general wavelength shift that tends to align the absorption features may bring out similari- 
ties both to the eye and in the DIFF2 test (i.e., by minimizing the DIFF2 function value). 



Branch et al.l (l2006al ) showed that general shifts were useful in bringing out similarities. We 
expect that the DIFF2 test shift will usually minimize the DIFF2 function value by tending 
to align blueshifted absorption features in photospheric phase spectra. 

In the nebular phase, the emission lines tend to be symmetric about the rest-frame 
line-center wavelength. If this is the case, the minimizing shift of the DIFF2 test is likely to 
be small or zero and DIFF2 will tend to reduce to DIFF1. However, if the emission peaks of 
one or both locally-normalized spectra are offset from the line-center wavelengths because 
of some asymmetry in the ejecta, then the DIFF2 test will compensate for that and again 
bring out similarities that could be missed by the DIFF1 test. 

Because DIFF2's logarithmic wavelength shift will often compensate for Doppler shifts 
caused by differences in velocity structure, it is useful to define a corresponding Doppler 
shift velocity parameter D param characteristic of the difference in velocity scale between two 
locally-normalized spectra. We first preliminarily define w pa ram for a wavelength shift 5X 
from an initial wavelength A using the 1st order Doppler formula (eq. (JSJ) in § 2.2, mutatis 
mutandis): thus we have 

oA ^param /1 o\ 

~r = > Uo) 

A c 

where we define positive velocities as giving redshifts. Say the corresponding logarithmic 
wavelength shift is 5ln(A) and recall we only allow logarithmic wavelength shifts of up to 
0.05. The allowed logarithmic wavelength shifts are sufficiently small that we can then make 
the approximation 

<fln(A)«^. (19) 
A 

We now give the final definition for t> param by combining equations ([TBI and (fl9j) : 

Wparam = c5 hl(A) . (20) 
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4.3. An Example of Two Locally-Normalized Spectra Fitted by Minimizing 

the DIFF2 Function Value 

Figure 6 shows an example of two locally-normalized spectra with the dotted-line spec- 
trum given the logarithmic wavelength shift that minimizes the DIFF2 function value for 
the two spectra. The solid-line spectrum is the SN 1993J 1993 April 15 spectrum (which 
comes from about 2 days before UVOIR bolometric maximum light) shown i n Figure 2 



3.2) and the dotted-line spectrum is the SN 1987K 1987 August 9 spectrum (IFilippenko 



19881 ) (which comes from about 9 days after a red-optical maximum light). In this case, 
DIFF1 = 0.642 and DIFF2 = 0.612 and the Doppler shift velocity parameter (see definition 
eq. (1201) in § 4.2) corresponding to DIFF2's logarithmic wavelength shift has absolute value 
999 km s -1 . In this case, DIFF2's logarithmic wavelength shift is a blueshift when applied to 
the SN 1987K spectrum: as noted above we have given this shift to the SN 1987K spectrum 
in Figure 6. The change from DIFF1 to DIFF2 and the Doppler shift velocity parameter 
are not large in this case and, in fact, the eye cannot detect any overall improvement in fit 
in going from a plot without the shift to the plot with it (i.e., Fig. 6). The small change 
between the DIFF1 and DIFF2 and the small shift is to be expected. First, the spectra 
both come from Type lib supernovae. Second, out of our current and preliminary sample 
of hydrogen-deficient spectra (see § 5), the SN 1987K 1987 August 9 spectrum is the closest 
match according to DIFF2 to the SN 1993J 1993 April 15 spectrum, except for other spectra 
from SN 1993J that come from phases close to 1993 April 15. According to DIFF1, it is only 
the second closest non-SN-1993J match: the SN 1987K 1987 August 7 spectrum is slightly 
closer with DIFF1 = 0.624. 

The phase (day —2) of the SN 1993J spectrum is relative to the well determined UVOIR 
bolometric maximum light (see §§ 3 and 3.2). The phase (day +9) of the SN 1987K spec- 
trum is relative to a red -optical maximum that is estimated to be 1987 July 31 with an 



uncertainty of ±4 days (jFilippenkd 119881 ). Supernova SN 1993 J was a particularly well- 
observed supernova, and hence our precise knowledge of its phase. Supernova SN 1987K was 
only moderately well-observed, and hence our uncertainty about its phase and also the light 
curve from which the phase was determined. The data available for SN 1987K is typical for 
the supernovae that must dominate any current sample of supernovae for statistical spectral 
analysis. 

The identifications in Figure 6 are the same as in Figure 2 (§ 3.2), except for the residual 



telluric absorption lines in the SN 1987K spectrum (jFilippenkd Il988l ) . The region from the 
blue side of the telluric lines to the red end of the SN 1987K spectrum were excluded from the 
evaluation of DIFF1/2: i.e., this region was not part of the good overlap region in logarithmic 
wavelength for the spectrum pair (§§ 4.1 and 4.2). 
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4.4. Comparing DIFF1 and DIFF2 



Both DIFF1 and DIFF2 measure the physical similarity of line formation, and thus of 
supernovae. DIFF1 puts more weight on overall physical similarity of the supernovae than 
DIFF2 because the DIFF2 shift in logarithmic wavelength compensates somewhat for some 
physical distinction in the supernovae. A spectrum pair with very small DIFF1 value may 
well come from nearly identical supernovae at nearly the same phase. On the other hand, a 
spectrum pair with a very small DIFF2 value, but not a small DIFF1 value, may well come 
from supernovae that are in some ways physically alike, but are physically distinct in some 
other ways. 

The most obvious physical distinction that the logarithmic wavelength shift in the 
DIFF2 test can compensate for is in velocity structure: the distinction is caused either by an 
intrinsic difference or a phase difference. Another possible distinction is in viewing direction 
for supernovae that are asymmetric. Note that asymmetric supernovae even if identical would 
have different spectral evolution depending on viewing angle. There is, in fact, considerable 
evidence from su pernova polarimetry and spectropolarimetry f or asymmetry in core-collaps e 



supernovae (e.g., iWheeler fc Benettil l2000l ; IWang et al.l l200ll ; iLeonard fc Filippenkd 120051 ) . 
One important kind of asymmetry can be parameterized by a major-to-minor axis ratio for 
characteristic orthogonal axes perpendicular to the line of sight. Spectro polarimetric obser - 
vations suggest that this axis ratio can vary from 1 to 2 or more (e.g., IWang et al.ll200ll ). 
For example, the Type lib supernovae SN 1993J and S N 1996cb spectrop olarimetric data 
suggest an axis ratio of > 1.4 for both these supernovae (jWang et al.l 1200 ll ). There are data 
that suggest, but with co nsiderable uncertainty, that asymmetry increas es with decreasing 
hydrogen envelope mass (I Wang et al.l l200ll ; ILeonard &: Filippenkd 120051 ). Another kind of 
asymmetry in core-collapse supernovae, that of a jet (or bipolar jets) emerging from the main 
part of the ejecta, is suggested by the spect ropolarimetry of Type Ic hypernova SN 2002ap 
( jKawabata et al.l 120021 ; ILeonard et al.l 120021 ) . Type la supernovae also show polarization i n 
some cases: the polarization may arise from clumps in the ejecta (e.g.. ILeonard et al.ll2005l ). 



Because of the use of local normalization, both DIFF1 and DIFF2 tests are useful for 
studying supernovae which frequently have uncertain continua. 



5. A PRELIMINARY STATISTICAL ANALYSIS OF 
HYDROGEN-DEFICIENT CORE-COLLAPSE SUPERNOVA SPECTRA 



We are beginning a project of comparative analysis of the spectra of HDCC supernovae 
(i.e., hydrogen-deficient core-collapse supernovae of Types lib, lb, Ic, and Ic hypernovae: see 
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[IJ that will include a statistical analysis of the spectra as well a spectral modeli ng; using 



the parameterized spectrum-synthesis code SYNOW (e.g.. iBranch et al.ll2003l . 120051 . and ref- 
erences therein). We note that Types lib, lb, and Ic seem to be a sequence of decreasing 
hydrogen and helium in the supernova ejecta. Type lb supernovae have conspicuous helium 
lines, but not conspicuous hydrogen lines. Type lib supernovae have conspicuous, but weak 
hydrogen lines and resemble Type lb supernovae. Type Ic supernovae lack conspicuous hy- 
drogen lines and conspicuous helium lines. As mentioned in the § HJ Type Ic hypernovae 
are particularly energetic Type Ic supern ovae that in some ca ses are associated with GRBs: 
SN 1998bw with GRB 980425 (e.g., balama et alJll998h . All the core-collapse super- 



e.g 



novae lack the strong Si II A6355 line that cha racterizes Type la supernovae. The spectral 
classification of supernovae is well reviewed by iFilippenkol (119971 ). but note that this paper 
precedes the discovery of Type Ic hypernovae. 

We need to make clear that the HDCC supernova types are not completely distinct and 
the determination of whether particular lines are conspicuous or not is sometimes ambiguous. 
For exampl e, relatively inconsp icuous hydrogen lines do seem to be present in Type lb super- 



novae (e.g., IBranch et al.ll2002l ). and so drawing the line between Type lib's and Type Ib's 
cannot always be clear. Likewise there is sometimes difficulty in distinguishing Type Ib's 
from Type Ic's and some supernovae that cannot be clearly set in e ither category ar e often 
classified as Type Ibc supernovae. On the other hand, the work of iMatheson et al.l (120011 ) 
shows that the distinction between Type Ib's and Type Ic's is strong enough to continue 
requiring the two types rather than a single merged type. 

To undertake the HDCC supernova analysis project, we are assembling HDCC supernova 
spectra in html files in supernova directories at SUSPEND (see § 1) under the heading 
Supernovae by Epoch. The spectrum files have header information about the spectra and 
the supernovae they come from. Plots of the spectra in original form and locally-normalized 
form are also given in the headers. The original spectra follow the header in two-column 
format: the columns being wavelength and flux. The locally-normalized spectra are given in 
the dif subdirectories of the supernova directories. A list of the spectra are currently in html 
files (including those spectra from non-HDCC supernovae) can be found under the heading 
Lists at SUSPEND. 

At present, we are at a very early stage in the project and have not assembled many 
of the spectra available for HDCC supernovae nor determined the best way to proceed with 
our statistical analysis. However, to illustrate the use of DIFF1 and DIFF2, we present some 
preliminary results. 



-33- 



5.1. Table 1 and Tables at SUSPEND 



Table [Hgives some statistics on the HDCC supernovae and spectra we have assembled so 
far in SUSPEND (see § 1). As one can see from the table, there are currently 17 supernovae 
and 168 spectra in the sa mple. The supernova e include the prototype Type lb supernovae 
SN 1983N and SN 1984 L Jrlarkness et al.lll987h and the well-observed Type Ic' s SN 1987M 
(IFilippenko et al.lll990l ) and pr ototypes. SN 19941 (e.g.. IClocchiatti et al.lll996l ). The pecu- 
liar Type Ibc SN 2005bf (e.g., iFolatelli et all hood : IParrent et alihoOTf ) is included as well: 
this supernova had two observed maxima, had a rise time of about 40 days to the second 
and main optical maximum (typical Type Ib/c rise times are of order 20 days), and was 
unusually luminous. (Supernovae also have a near-explosion maximum light that is almost 
never observed, and so usually unmentioned. This is not the maximum light we refer to here 
as the first maximum light of SN 2005bf.) SN 2005bf also made a transition from Type Ic 
to Type lb behavior in the pe riod from about day —30 to about day +20 relative to the 



(main) optical maximum (e.g 
is classified clS db Type Ibc. Only two Type II 
already discussed prototype S N 1987K (§ 4 .3; 
SN 1993J (§ 3 and S 4. 3: e.g.. Ijefferv et al. 



Folatelli et al.ll2006l ): because of this transition, SN 2005bf 



3 supernovae hay e been included so far: the 



Clocchiatti et al. 1995: Houck fc Fransson 



pernovae are included: SN 



1998: 



Iwamoto et al. 



Leonard et al 



200^ 



199 



997ef (e.g., 



Patat et a 



observed by iMatheson et al.l (120011 ) and others 



Fi 



1994 



1996: 



ippenkolll988r) and the very well-observed 



Richmond et al. 1994: Baron et al. 1995 



Fransson et al.l 120051 ) . Thre e Type Ic hy- 



Mazzali et al.ll2000h . SN 1998 bw fe.g-.lGalama et al 



200 if ), and SN 2002ap (e.g., Kawabata et al. 



2002 



2); Foley et al. 2003 ). Also included are less famous HDCC supernovae 



In the calculations described below, we included only one spectrum for each phase period 
(which we usually set to being 1 day) from each supernova. The other spectra for that phase 
period are redundant since they should be nearly identical to the included spectrum aside 
from observational error. The spectrum we included is the one we deem to be best: this is 
usually because of broader wavelength coverage. 

In an automated fashion, we have calculated DIFF1 and DIFF2 for all supernova spec- 
trum pairs in the sample: these spectrum pair DIFF1/2 values can be found tabulated at 
SUSPEND under the heading Lists. For each spectrum, there is a table listing all other 
spectra in order of increasing DIFF1/2 values relative to it. These tables are updated as 
more supernovae and spectra are added to the sample. We have excluded from the analysis 
below and in § 5.2 those DIFF1 and DIFF2 values deemed invalid by the rules discussed in 
55 4.1 and 4.2. 



As well as ordering spectra in order of likeness using DIFF1/2 values, we want some way 
to order supernovae in order of likeness. As a preliminary method, for any two supernovae 
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we calculate all (valid) DIFF1/2 values between spectrum pairs with one of the pair coming 
from one supernova and the other from the other supernova: i.e., we calculate the non-self 
spectrum pair DIFF1/2 values. The smallest of these DIFF1/2 values we call the supernova 
pair DIFF1/2 value. 

The rationale for the definition for supernova pair DIFF1/2 value is as follows. The 
spectral phase coverage for all supernovae is incomplete (although for some the coverage is 
very good) and rarely are both supernovae covered at a common phase to within a day: often 
coverage at even a nearly common phase is lacking. Also phase measured from a time zero 
set at maximum light (which is how we define phase following a common convention) is itself 
uncertain. Ideally, one would like to use an accurate UVOIR bolometric maximum light as 
time zero for supernova phase. We do that when an accurate UVOIR bolometric maximum 
light is available (which is rarely). Usually, however, we must use a substitute which is often 
some optical broad-band maximum light. These broad-band maxima usually happen within 
a few days of UVOIR bolometric maximum light and have some uncertainty. In some cases, 
all we can use as a substitute is what can be loosely called optical maximum light which is a 
fiducial time within a few days of when exact UVOIR bolometric maximum light probably 
occurs. And in some of these cases, the estimate of optical maximum light is just the day 
of discovery: supernovae tend to be discovered near or at optical maximum light since that 
is when they are most readily discoverable because, of course, they are brightest or nearly 
brightest in all optical bands at optical maximum light. (Throughout the rest of this section 
we usually do not specify the kind of maximum light we use to set the phase zero point. We 
usually just refer to all kinds of maximum light as maximum light without qualification.) 
By choosing the smallest DIFF1/2 value out of all the non-self spectrum pair values for the 
supernova pair DIFF1/2 value, we hope to have partially avoided the problems of incomplete 
spectral phase coverage and inaccurate phase measured from UVOIR bolometric maximum 
light. We assume that the intrinsic likeness of each supernova pair is well measured by 
their closest observed spectral approach to each other: i.e., what we call the supernova pair 
DIFF1/2 value. Obviously, because of the incomplete phase coverage, the closest observed 
spectrum approach may only give a lower limit on likeness. For two supernovae that are 
actually much alike with at least one having good phase coverage, the closest spectrum 
approach is probably for the spectrum pair that has the least relative phase difference where 
in this case phase is measured relative to explosion. (As supernovae evolve from time zero 
at explosion, their relative rate of spectral evolution tends to be more constant than their 
absolute rate of spectral evolution.) For two supernovae that are not much alike or that both 
have poor phase coverage, the closest spectrum approach may not occur for the spectrum 
pair that has the least or even small relative phase difference. 

There are, in fact, cases where two supernovae have a relatively small supernova pair 
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DIFF1/2 value and yet are in some respects quite different. For example, the day +4 (relative 
to B maximum) spectrum of Type Ibc SN 1999ex ( Hamuy et al.ll2002l ; IStritzinger et al.ll2002l ) 



is very similar to the day —20 (relative to UVOIR bolometric maximum light ) spectrum of 



Type Ibc SN 2005bf flFolatelli et al.ll2006l : IParrent et al.ll2007h . Figure 10 of IParrent et al. 



(120071 ) shows overlapping loca lly- normalized versio ns of these spectra: to the eye their fit is 



very close. (See also Fig. 5 of iFolatelli et al.ll2006l ) The supernova pair DIFF1 and DIFF2 



values are, respectively, 0.447 and 0.442: these are quite small DIFF1/2 values for spectra 
from different supernovae as our experience shows. Although SN 1999ex and SN 2005bf are 
both Type Ibc supernovae, they are quite different in so me respects: e.g., rise tim e to optical 
maximum: probably of order 20 days for SN 1999ex (IRichardson et al.l 120061 ) and about 
40 days to the second and main optical maximum, as mentioned above, for SN 2005bf. Also 
the two closely matching spectra come from very different phases which would not be true if 
the supernovae were really very much alike: the most closely matching spectra would then 
tend to come from nearly the same phase. Thus, it is clear that our assumption of intrinsic 
supernova likeness being well measured by the supernova DIFF1/2 pair is not always valid. 
However, since our analysis here is preliminary, we leave the search for a better measure of 
supernova likeness to future work sine die. 

The tables giving the spectrum pair DIFF1/2 values at SUSPEND (which altogether are 
quite lengthy) are preceded by tables of supernova pair DIFF1/2 values. For each supernova 
there is a table listing all other supernovae in order of increasing supernova pair DIFF1/2 
values relative it. 



5.2. Table 2 and the Standard HDCC Supernova Types 



Since this paper is just a beginning in our project of statistical analysis of HDCC super- 
nova spectra, we will not fully discuss the significance of either the tabulated supernova or 
spectrum pair DIFF1/2 values at SUSPEND (see § 5.1). But as a preliminary investigation, 
we consider if the classification of HDCC supernovae into Type lib's, Type Ib's, Type Ic's, 
and Type Ic hypernovae is discoverable from supernova pair DIFF1/2 values. (Supernovae 
classified as Type Ibc are g rouped with the Type Ib's. There are onl y two of these in the 
current sample: SN 1999ex rtHamuv et alJliooilstritzinger et alJbooj ) and SN 2005bf (e.g., 
Folatelli et al.ll2006l ; IParrent et al.ll2007l ).) A priori, it is not clear that the HDCC supernova 
types will be picked out by an analysis with DIFF1/2. The Types lib, lb, and Ic, and, in 
part, the Type Ic hypernova are characterized by particular lines. DIFF1/2 measures overall 
differences between spectra. 



For the preliminary investigation, we have calculated the mean supernova pair DIFF1/2 
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values between supernovae of different types (the cross-type means) and the mean supernova 
pair DIFF1/2 values between supernovae of the same type (the self-type means). We have 
also calculated the estimated standard deviations of the distributions of the supernova pair 
DIFF1/2 values about the type means using the ord inary standard d eviation formula with 



the correction for using the mean of the sample (e.g.. lBevingtonlll969l . p. 19). The standard 
deviation for the Type lib self-type mean is assigned a zero value since there are only two 
Type lib's in the current sample, and thus only one supernova pair DIFF1/2 value: the 
standard deviation of the sample is zero; that of the distribution is unknown. 

Since HDCC supernovae tend to look increasingly alike as time from explosion (i.e., 
phase) increases, we have repeated the DIFF1/2 calculations including only supernova pair 
DIFF1/2 values obtained from spectra that are both from phases less than 10 days past 
maximum light. The supernova pair DIFF1/2 type means obtained in this case are phase- 
restricted. We were interested in seeing if phase-restricted approach would better distinguish 
HDCC supernovae into the standard types than the phase unrestricted approach. Note that 
phase-unrestricted supernova pair DIFF1/2 values will, of course, always be smaller than or 
equal to their counterpart phase-restricted supernova pair DIFF1/2 values. (The last remark 
assumes that counterpart phase-restricted supernova pair DIFF1/2 values exist. They may 
not if one or both of the supernovae in the pair lack any spectra for the period allowed by 
the phase restriction.) 

The results of all the calculations are given in Table 2 in the form of a type DIFF1/2 
correlation matrix consisting of self-type means (the diagonal elements) and cross-type 
means (the off-diagonal elements): phase-unrestricted DIFF1 values are in Table 2(a), phase- 
unrestricted DIFF2 values are Table 2(b), phase- restricted DIFF1 values are in Table 2(c), 
and phase-restricted DIFF2 values are in Table 2(d). The matrix is, of course, symmetric. 

Several remarks can be made about the results in Table 2. First, none of the mean 
values are very small: the smallest is 0.455 which is the Type Ic hypernova DIFF2 self-type 
mean in Table 2(b). An examination of the table of the phase-unrestricted supernova pair 
DIFF1/2 values at SUSPEND shows that none are less 0.4: the smallest DIFF1 is 0.418 for 
the pair of Type lb SN 1988L and Type Ic SN 1990B; the smallest DIFF2 is 0.409 for the 
same pair of supernovae and obtained using the same spectrum pair. The fact that these 
smallest DIFF1/2 values are for supernovae of different assigned types shows that supernova 
pair DIFF1/2 values are not perfect in discriminating Type lb and Type Ic supernovae. 
Some of the spectrum pair DIFF1/2 values are significantly smaller (some below 0.3 going 
down to minimum values of 0.244 for DIFF1 (see § 4.1) and 0.229 for DIFF2 (see § 4.2)), 
but those are for spectra from the same supernova usually taken close together in phase. 
The upper limit on the phase-unrestricted supernova pair DIFF1 values is 0.886 and on the 
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phase-unrestricted supernova pair DIFF2 values is 0.809. These upper limits were both for 
the supernova pair of Type Ic SN 1990aa and Type lb SN 1999di, but not using the same 
spectrum pair. The supernovae of this pair are not much alike in the phases for which spectra 
are available. 

Second, a mean DIFF2 value is always less than or equal to the corresponding mean 
DIFF1 value as it must be according to our formulae (§ 4.1, eq. (12) and § 4.2, eq. (17)), but 
always by less than the standard deviation of the mean DIFF1 value. We have to conclude 
that DIFF2 is not finding great similarities that are being totally missed by DIFF1. 

Third, the phase-restricted mean DIFF1/2 values (in Table 2(c) and Table 2(d)) are 
all larger than the corresponding phase-unrestricted mean DIFF1/2 values (in Table 2(a) 
and (b)). This is understandable for two causes. First, HDCC supernovae are qualitatively 
understood to be more diverse at earlier times, and thus the phase-restricted DIFF1/2 values 
should tend to be larger. Second, the incompleteness of the phase coverage of a supernova 
may be greatly increased by phase restriction. Phase restriction reduces the chances of finding 
supernova pair DIFF1/2 values for the aforementioned supernova with another supernova 
with pair spectra from nearly the same phase when the spectra would tend to be most alike 
if the pair supernovae were alike. The effects of the two causes could be sorted out in a 
more detailed analysis. It is possible for the mean DIFF1/2 values to be reduced by phase 
restriction. This would happen if phase restriction eliminated a supernova altogether from 
those going into calculation of the means and that supernova was sufficiently remote from 
other supernovae to increase phase- unrestricted mean DIFF1/2 values when included. This 
actually happened in an earlier calculation for the Type Ic self-type means. The phase- 
restriction eliminates Type Ic SN 1990aa for which no spectra from before 10 days after 
maximum light are currently available to us. SN 1990aa may be something of an outlier 
among Type Ic's, and so its inclusion seems to increase the self- type mean DIFF1/2 values. 
When the sample of Type Ic's was changed this reduction on phase restriction disappeared. 

Fourth, the self-type means in any row are the smallest means in that row in all cases. 
This shows that a DIFF1/2 analysis can probably identify the standard HDCC supernova 
types on average. If one is given spectra from a new HDCC supernova of unknown type and 
uses those spectra to find mean supernova pair DIFF1/2 values with samples of the standard 
types, the smallest mean supernova pair DIFF1/2 value will probably be with the sample of 
the type to which the new HDCC supernova belongs. 

In the rows for Type lib's and Type Ic hypernovae, the self-type means are smaller 
than the respective row cross- type means by > 1.5cr larger in all cases and by > 2<Ti arger in all 
but two cases. (The o"i argC r is the larger of the standard deviations for the two means being 
compared.) The two cases where difference is less than 2o"i argcr are for the DIFF1 and DIFF2 
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phase-unrestricted Type Ic hypernovae cross-type with the Type Ic's. Overall, Type lib and 
Type Ic hypernova are found to be quite distinct from other types. 

On the other hand, in the phase-unrestricted rows for Type Ib/Ibc's and phase- unrestricted 
and phase-restricted rows for Type Ic's the self-type means are smaller than the respective 
row cross- type means by < o-i arger . For the phase- restricted cases, the Type Ib/Ibc self- type 
means in are some cases separated by more than o"i arg cr from the Type Ib/Ibc row cross- 
type means. Overall, one must conclude that Type Ib/Ibc and Type Ic are not as distinct 
from other types as are Type lib and Type Ic hypernova. This means that Type Ib/Ibc 
and Type Ic are more heterogeneous and overlap more in behavior with other types than 
Type lib and Type Ic hypernova. 

Recall that we introduced phase restriction to see if it would better distinguish HDCC 
supernovae into the standard types than the phase unrestricted approach. From the results 
discussed in the last two paragraphs, there is some improvement in distinguishing the types. 
However, this improvement is rather modest for the particular phase restriction we used: 
i.e., restricting spectra used to those from before 10 days after maximum light. Maybe other 
phase restrictions would give greater improvement. This could be investigated in future 
work. 

To finish our preliminary analysis, we note that the Type lib's and Type Ic hypernovae 
included in the sample are all relatively well-observed supernovae, and so they suffer relatively 
little from lack of phase coverage. Some of the Type Ib/Ibc's and Type Ic's in the sample have 
poor or very poor phase coverage. (See the list of supernova spectra in html format under 
the heading Lists at SUSPEND for the spectra currently available for individual supernovae 
and the spectrum phases.) The supernova pair DIFF1/2 values for supernovae where both 
supernovae have poor phase coverage could be large even if the supernovae were very similar 
simply because the procedure for determining the supernova pair DIFF1/2 values may not 
be able to find spectra from nearly the same phase. In future work, we will try to remedy 
this deficiency in the procedure perhaps by giving lower weight in calculating mean DIFF1/2 
values to supernova pair DIFF1/2 values from supernova pairs where both members of the 
pair have poor phase coverage or by using some quite different definition of supernova pair 
DIFF1/2 value. 

6. CONCLUSIONS AND DISCUSSION 

We have developed two tests DIFF1 and DIFF2 for measuring goodness-of-fit between 
two supernova spectra (see § 4.1, eq. (12) and § 4.2, eq. (17) for the formulae). The tests rely 
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on local normalization (§ 2) which eliminates the uncertainty in the spectrum continuum, 
except for contamination from extraneous sources (§ 1). Local normalization also eliminates 
real information stored in continuum shape. However, a basic premise of this paper is 
that line pattern is a much better signature of intrinsic supernova behavior than continuum 
shape, and so eliminating continuum shape information is not too important for spectrum 
comparisons (§1). 

We have presented some examples of locally-normalized spectra for SN 1993 J and given 
some analysis of the spectra (§ 3). The UV parts of two of the SN 1993 J spectra are hitherto 
unpublished HST spectra. In § 3.5, we have shown that local normalization is sufficiently 
continuum-independent to lead to valid results in comparing locally normalized spectra. One 
must, however, apply the same local normalization procedure to all original spectra in order 
to obtain locally-normalized spectra that can be validly compared. 

As an example of the use of DIFF1/2, we have used them in a preliminary statistical 
analysis of the spectra of HDCC (hydrogen-deficient core-collapse) supernovae (§ 5). The 
analysis shows that standard HDCC supernova types (lib, Ib/Ibc, Ic, and Ic hypernova) do 
form distinct groups when compared using DIFF1/2 (§ 5.2). This analysis is preliminary 
since many available HDCC supernovae and HDCC supernova spectra are not included in 
the analysis sample. Also many improvements in our statistical procedure are possible. 

Tables of all spectrum and supernova pair DIFF1/2 values for our sample of HDCC 
supernovae are available at the SUSPEND database (see the footnote to the abstract for the 
URL) under the heading Lists. The spectra we have used are also online at SUSPEND in 
two-column format under the heading Supernovae by Epoch in supernova directories in html 
files along with figures of the original and locally-normalized spectra. The locally- normalized 
versions of the spectra themselves can be found in the dif subdirectories of the supernova 
directories. 

Although hu ndreds of supernovae are now being discovered per year (e.g., 367 in 2005 



and 526 in 2006 ([Central Bureau for Astronomical Telegrams! 120071 ) which gives discovery 



rates of more than 1 supernova per day), most of these are remote and are relatively poorly 
observed. Their main use is as cosmological distance indicators and for the determination 
of supernova rates. New well-observed supernovae accumulate slowly with only a few per 
year. These new well-observed ones and the past well-observed ones are only relatively 
well-observed in most cases. Phases are missing, calibrations imperfect, and frequently the 
reddening correction is very uncertain. Thus, for the foreseeable future, statistical analyses 
of supernova spectra will have to rely on heterogeneous data sets for relatively few well- 
observed supernovae. We believe that DIFF1 and DIFF2 (which eliminate the need for 
accurate continuum shape and that can, in the case of DIFF2, compensate somewhat for 
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varying phase, velocity structure, and asymmetry (§ 4.4)) will be useful tools in statistical 
analyses of available spectra. 

Another, and not-distinct, use for DIFF1/2 is to find for new and, perhaps, not-well- 
observed supernovae, well-observed supernovae that are their near-twins. Then insofar as 
the well-observed supernovae are understood, the new supernovae will be understood. 

DIFF1/2 should also be useful in synthetic spectrum analysis of supernova spectra 
because, again, it eliminates the uncertain continuum of observed supernovae. Certainly 
fitting the true continuum is one of the goals and guides in achieving a good synthetic 
spectrum fit to the observations. But when the fitting exercise gives a good fit to an incorrect 
continuum, then the fitting exercise becomes misleading. A good fit to the lines with highly 
realistic radiative transfer and a realistic hydrodynamic model should yield a good continuum 
and that should allow one to correct the observed continuum for errors in dereddening and 
calibration. 

The future of synthetic spectrum modeling of supernovae may well be the calculation 
of time-dependent radiative transfer for a large collection of realistic hydrodynamic models 
that span many possible explosion outcomes. As time passes more models will be included 
in the collection and the realism of the models and the radiative transfer will be improved. 
Eventually, almost all old and new supernovae could find a near-twin in the collection using 
some test like DIFF1/2 and will then be well understood. Thi s ideal situation is still far off, 



but there is work heading toward it (e.g., IWoosley et al.l 120071 . for Type la supernovae) 



In future work, we plan to use DIFF1/2 in fuller statistical analyses of HDCC supernova 
spectra than that presented in this paper. These analyses will also include spectral modeling 



using the parameterized spectrum-synthesis code SYNOW (e.g., iBranch et al.l 12003k 12005 
and references therein). Use of DIFF1/2 for other supernova types is also envisaged for 
future work. 
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referee for his/her comments. Support for this work has been provided by NASA grant 
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Fig. 1. — The spectrum of Type lib SN 1993J in the fx and locally-normalized representa- 
tions from 1993 March 31 which is about 16 days before UVOIR bolometric maximum light 
and about 4 days after explosion. The locally-normalized spectra in this and in other figures 
are obviously the ones with continuum level of about 1. 
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Fig. 2. — The spectrum of Type lib SN 1993J in the fx and locally-normalized represen- 
tations from 1993 April 15 which is about 2 days before UVOIR bolometric maximum light 
and 18 days after explosion. From 3240 A blueward, the spectrum is an HST spectrum. 
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Fig. 3. — The spectrum of Type lib SN 1993J in the f\ and locally-normalized represen- 
tations from 1993 September 15.5 (an averaged date) which is about 151 days after UVOIR 
bolometric maximum light and 171 days after explosion. From 4253.6 A blueward, the spec- 
trum is an HST spectrum (from 1993 September 17). 
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Fig. 4. — The spectrum of Type lib SN 1993J in the f\ and locally-normalized represen- 
tations from 1993 November 15.5 (an averaged date) which is about 212 days after UVOIR 
bolometric maximum light and 232 days after explosion. From 4200 A blueward, the spec- 
trum is an HST spectrum (from 1993 November 14). The original HST spectrum blueward 
of ~ 2425 A seems to decline in a manner too steep to be physically real and we do not trust 
it nor the locally-normalized spectrum there. 
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Fig. 5.— The evolution of the blueshifted Mg II A2797.9 emission line in the locally- 
normalized spectra of SN 1993J. The phases given in parentheses in the figure are relative to 
the UVOIR bolometric maximum light on 1993 April 17. The horizontal axis is logarithmic, 
but this is hard to notice given the small wavelength range covered. 
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Fig. 6.— The locally-normalized spectra of Type lib supernovae SN 1993J from 1993 
April 15 (about 2 days before UVOIR bolometric maximum light) and SN 1987K from 1987 
August 9 (about 9 days after a red-optical maximum light). The SN 1987K spectrum has 
been blueshifted using Doppler shift velocity parameter absolute value 999 km s -1 to mini- 
mize the DIFF2 function value for the spectrum pair. For the spectrum pair, DIFF1 = 0.642 
and DIFF2 = 0.612. The spectra were both initially corrected for host galaxy heliocen- 
tric velocity: from Leda, the mean host gal axy heliocentric ve locities are —39 ± 2kms _1 
(SN 1993J) and 799±2kms" 1 (SN 1987K) JPaturel et alil2003h . 
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Table 1. Current Sample in SUSPEND of HDCC Supernovae and Spectra 



Type 


Number of Supernovae 


Number of Spectra 


All HDCC Types 


IV 


168 


Type lib 


2 


39 


Type Ib/Ibc 


8 


45 


Type Ic 


1 


34 


Type Ic hyp 


3 


50 



a Type Ibc supernovae are those supernovae which are not clearly 
distinguishable into Type Ib's or Type Ic's. We group Type Ib's and 



SN 1999ex jHamuv et al. 


|2002|; IStritzineer et alj|2002l) and SN 2005bf 


(e.g.,lFolatelli et al .1120061; 


Parrent et al.l 20071). Tvdc Ic hvp is short for 



Type Ic hypernovae. 
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Table 2. HDCC Supernova Type DIFF1/2 Correlation Matrix 



TypeVType 


Type lib 




Type Ib/Ibc a 




Type Ic 




Tvpe Ic hvpk 






Mean DIFF1/2 


St. Dev. 


Mean DIFF1/2 


St.Dev. 


Mean DIFF1/2 


St. Dev. 


Mean DIFF1/2 


St.Dev. 




(a) DIFF1: 


no phase restriction 












Type lib 


0.507 


0.000 


0.674 


0.077 


0.662 


0.068 


0.632 


0.052 


Type Ib/Ibc 


0.674 


0.077 


0.615 


0.097 


0.646 


0.113 


0.712 


0.081 


Type Ic 


0.662 


0.068 


0.646 


0.113 


0.578 


0.118 


0.614 


0.085 


Type Ic hyp 


0.632 


0.052 


0.712 


0.081 


0.614 


0.085 


0.499 


0.068 




(b) DIFF2: 


no phase restriction 












Type lib 


0.507 


0.000 


0.644 


0.071 


0.618 


0.038 


0.601 


0.042 


Type Ib/Ibc 


0.644 


0.071 


0.573 


0.094 


0.610 


0.105 


0.670 


0.056 


Type Ic 


0.618 


0.038 


0.610 


0.105 


0.547 


0.120 


0.579 


0.071 


Type Ic hyp 


0.601 


0.042 


0.670 


0.056 


0.579 


0.071 


0.455 


0.020 


(c) DIFF1: restricted to pre-10-days past maximum light 


Type lib 


0.507 


0.000 


0.730 


0.064 


0.699 


0.066 


0.735 


0.086 


Type Ib/Ibc 


0.730 


0.064 


0.633 


0.085 


0.654 


0.082 


0.855 


0.128 


Type Ic 


0.699 


0.066 


0.654 


0.082 


0.597 


0.160 


0.754 


0.135 


Type Ic hyp 


0.735 


0.086 


0.855 


0.128 


0.754 


0.135 


0.504 


0.072 


(d) DIFF2: restricted to pre-10-days past maximum lif 












Type lib 


0.507 


0.000 


0.696 


0.053 


0.647 


0.038 


0.681 


0.069 


Type Ib/Ibc 


0.696 


0.053 


0.593 


0.076 


0.623 


0.071 


0.763 


0.074 


Type Ic 


0.647 


0.038 


0.623 


0.071 


0.557 


0.130 


0.652 


0.083 


Type Ic hyp 


0.681 


0.069 


0.763 


0.074 


0.652 


0.083 


0.485 


0.064 



a Type Ibc supernovae are those supernovae which are not clearly distinguishable into Ty pe Ib's or Type Ic's. We group the 
Type Ib's and Type Ibc's together. There are only two Type Ibc's in the sample: SN 1999ex l lHamuv et al.l l2002: Stritzingc r et all 
l2002h and SN 2005bf (e.g.. iFolatelli et al.|[200ri IParrent et al.ll2007l) . 

b Type Ic hyp is short for Type Ic hypernovae. 

Note. — This table gives the mean DIFF1/2 values between HDCC supernova types in a correlation matrix format. The DIFF1/2 
value between two supernovae in this preliminary analysis is just defined to be the smallest DIFF1/2 value found out of all the pairs 
of spectra with one spectrum drawn from the one supernova and the other from the other supernova. With no phase restrictions 
the spectra are drawn from all phases. With phase restriction, the spectra are only drawn from the phases allowed by the given 
restriction. (The kind of maximum light used in setting the phase zero point is chosen various ways as discussed in § 5.1.) The 
matrix is symmetric, of course. The diagonal elements are the mean DIFF1/2 values for pairs with supernovae of the same type in 
the pairs: we call these means the self- type means. The off-diagonal elements are the mean DIFF1/2 values for pairs with supernovae 
of different types in the pairs: we call these means the cross-type means. More digits are shown than are significant to allow for 
numerical consistency checks and to give the table a consistent appearance. 



